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1.  General  Remarks. 

During  the  oeriod  November  1,  1981  to  October  31,  1982,  the 
Principal  Investigator,  in  cooperation  with  several  research 
assistants,  carried  out  a  urogram  of  mathematical  research  in 
the  general  area  of  control  theory  of  partial  differential  eaua- 
tions.  The  program  involved  two  distinct  phases?  an  effort  aimed 
specifically  at  the  development  and  imorovement  of  control  strsfc- 
egies  in  connection  with  the  wing  flutter  problem  and  a  more 
general  program  in  the  area  of  distributed  parameter  control  prob¬ 
lems  art  hyperbolic  type. 

This  work  resulted  in  two  scientific  papers  which  form  the 
greater  part  of  this  report..  The  first  of  these,  "Some  Remarks 
on  the  Current  Status  of  the  Control  Theory  of  Single  Soace  Dim¬ 
ension  Hyperbolic  Systems"  was  presented  at  the  NASA  JPL  Symposium 
on  Control  and  Stabilization  of  Large  Space  Structures,  Pasadena, 
CA,  July,  1982.  The  second,  "Admissible  Input  Elements  for  Systems 

in  Hilbert  Space  and  a  Carleson  Measure  Criterion",  by  L.  P.  Ho 

» 

and  the  Principal  Investigator,  is  a  paper  which  largely  resulted 
from  Dr.  Ho's  thesis  work,  also  supported  by  this  grant,  in  part. 

In  addition  to  Dr.  Ho,  who  is  now  with  the  University  of  Iowa, 
the  Principal  Investigator  was  assisted  by  R..  G.  Teglas,  H.  M. 


Baron,  and  R.  Rebarber.. 
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mot  funds  were  used  to  support  travel  by  the  Principal 
iHg* tor  and  one  Research  Assistant,  H.  U.  Baron. 

£«  Principal  Investigator  travelled  to  Pasadena,  California, 
end  and  present  a  paper  at  the  NASA  Jet  Prooulsion  Laboratory 
turn  on  Control  and  Stabilization  of  Large  Space  Structures. 
i®r  presented  concerned  the  current  status  of  the  control 
of  hyperbolic  partial  differential  equations  with  particular 
s  on  observers  and  canonical  structure. 

ie  Principal  Investigator  also  took  part  in  the  30th  Anni- 
’  meeting  of  the  Society  for  Industrial  and  Applied  Mathe- 
(SIAM)  in  Palo  Alto,  California,  July,  1982.  A  paper  out- 
the  treatment  of  control  problems  associated  with  infinite 
onal  linear  systems  by  means  of  methods  from  the  theory  of 
c  functions  was  presented  at  this  meeting  by  invitation  of 
;anizing  committee..  Ms.  Baron  also  attended  this  meeting 
sented  a  paper  on  control  canonical  forms  for  systems  gov- 
ly  various  types  of  partial  differential  equations. 


3 


3-  Technical  Appendix. 

This  appendix  consists  of  two-  papers  whose  preparation  was 

supported  in  part  by  the  grant.  These  papers  are: 

"Some  Remarks  on  the  Current  Status  of  the  Control  Theory 
of  Single  Space  Dimension  Hyperbolic  Systems'* 

and 

"Admissible  Input  Elements  for  Systems  in  Hilbert  Space 
and  a  Carleson  Measure  Criterion", 

the  latter  paper  being  jointly  authored  by  L.  F.  Ho  and  the 

Principal  Investigator. 


50 ME  REMARKS  ON  THE  CURRENT  STATUS  OF  THE  CONTROL  THEORY 
OF  SINGLE  SPACE  DIMENSION  HYPERBOLIC  SYSTEMS 


D.  L.  Russell’ 

University  of  Wisconsin  Mathematics  Department 

ABSTRACT 

We  review  various  aspects  of  the  control  theory  of  hyperbolic  systems, 
including  controllability,  stabilization,  control  canonical  form  theory,  etc.  To 
allow  a  unified  and  not  excessively  technical  treatment,  we  restrict  attention 
to  the  case  of  a  single  space  variable;  the  multi -dimensional  case  is  treated  in 
our  more  extensive  review  [  36]  .  The  paper  concludes  with  a  short  discussion 
of  the  newly  developed  procedure  of  canonical  augmentation. 

SOME  ASPECTS  OF  THE  CONTROL  THEORY  OF  THE  WAVE  EQUATION 

AND  RELATED  SYSTEMS 

The  systematic  study  of  control  systems  governed  by  partial  differential 
equations,  a  special,  but  exceptionally  important,  subcategory  of  distributed 
parameter  systems  began  in  the  early  1960's  with  the  work  of  the  Soviet 
scientists  A.  G.  Butkovskii  [  3  ]  ,  [4],  Yu.  V.  Egorov  [11]  and  others.  These 
works  were  primarily  concerned  with  the  extension  of  Pontyagin' s  Maximum 
Principle  [26]  to  certain  classes  of  processes  which  could  not  be  satisfactorily 
modelled  by  finite  dimensional  mathematical  systems.  Controllability  questions 
were  raised  but  were  usually  subsidiary  to  questions  of  optimality.  One  of  the 
first  systematic  controllability  studies,  in  connection  with  the  heat  equation, 
was  presented  by  Gal'chuk  in  [  14]  .  One  of  the  most  important  of  the  early 
American  contributions  to  the  subject  was  the  1963  thesis  of  Fattorini  [13]  , 
which  also  treated  parabolic  systems  and  was  one  of  the  first  works  to 
recognize  the  strong  relationship  between  distributed  parameter  control  studies 


*  Supported  in  part  by  the  Air  Force  Office  of  Scientific  Research  under  Grant 
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and  classical  results  in  analytic  function  theory. 

The  author*  s  own  interest  in  distributed  parameter  control  theory  arose 
out  of  consulting  experience  with  Honeywell,  Inc. ,  and  NASA,  starting  around 
1965  or  1966.  In  developing  the  Saturn  launch  vehicle  for  the  Apollo  program, 

NASA  has  encountered  the  problem  of  transverse  vibrations  of  the  booster 
structure  and  Interaction  of  those  vibrations  with  liquid  sloshing  modes  in  the 
immense  Saturn  fuel  tanks.  While  the  eventual  treatment  of  that  problem  was 
based  on  finite  modal  approximations,  the  problem  stimulated  a  great  deal  of 
research  aimed  at  an  understanding  of  the  control  of  vibrations  in  various 
distributed  parameter  settings. 

First  looking  at  this  problem,  under  Honeywell- NASA  auspices,  we 
thought  of  modelling  the  booster  structure  as  an  "Euler*1  beam,  the  displacement 
w  (x,  t),  which  we  may  take  to  be  scalar  here,  satisfying 

p(x)  +-4-  (EI(x)  ^  )  =  0  (1.1) 

at*  ax*  dxc 

along  with  appropriate  boundary  conditions  including  the  control  inputs,  at  the 
longitudinal  extremities  x  =  0.  x  =  L  .  We  got  nowhere  with  our  study  of 
this  problem  initially  because  the  equation  (1. 1)  is  not  particularly  well 
understood  from  the  mathematical  standpoint.  There  seemed  to  be  no  "handles" 
to  grasp,  ft  would  not  be  until  the  1969  thesis  of  Quinn  [27]  that  we  would 
understand  how  this  system  works  and  that  it  is,  in  fact,  controllable  in  a 
rather  strong  sense. 

We  knew  about  the  control  theory  of  ordinary  differential  equations  from 
various  papers  and  from  notes  and  lectures  which  would  later  be  incorporated 
into  the  1967  treatise  on  control  theory  by  Lee  and  Markus  [  20]  .  We  also 
knew  that  hyperbolic  partial  differential  equations  in  two  independent  variables 
reduce  to  ordinary  differential  equations  satisfied  along  the  characteristics,  ft 
was  natural,  therefore,  to  look  for  hyperbolic  models  which  might  fit  our  purpose. 
Such  was  provided  by  the  Timoshenko  beam  equations 

Iow^-kW(&-+,-£(EI(x,^,=  0  (1-2> 
r  8t 

pm  M  -  &«*»*•!£  -  ‘•'H =  °-  <*•*> 

dt 

which  may  be  viewed  as  two  coupled  wave  equations.  By  "wave  equations" 
here,  we  mean  the  equation 


-  b  <s « ■& )  ■ 


(1.4) 


All  coefficient  functions  shown  in  (1.  2),  (1.  3),  (1.  4)  are  positive  on  OsxsL. 

It  may  be  verified  that  (1. 2),  (1. 3)  and  (1. 4)  are  hyperbolic  in  the  sense 
described  in  [8],  [25],  e.g. .  Since  (1.4)  is  conceptually  simpler,  it  was 
studied  first,  with  accompanying  boundary  conditions 

z  (0,t )  =  0  (1.5) 


(1.6) 


the  latter  incorporating  the  control  force  u  (t)  . 

While  the  practical  goal  in  mind  was  appropriate  form  of  stabilization, 
we  knew  that  in  the  case  of  finite  dimensional  systems 

x  =  Ax  +  Bu 

an  affirmative  resolution  of  the  controllability  problem,  steering  from  a  given 
x(0)  =  Xq  to  a  given  x(T )  =  Xj,  Implied  the  property  of  stabilization;  hence 
we  felt  Justified  in  first  looking  at  the  state  to  state  controllability  problem  for 
(1. 4),  (1. 5),  (1.  6).  The  "energy"  form  for  (1. 4)  is 


(1. 4),  (1. 5),  (1.  6).  The  "energy"  form  for  (1. 4)  is 

t{z.  W  >  =  I  fQ  r M  <  TE  t »2  +  s<x>  ( IF  1 )  >2dx 

Given  initial  and  terminal  states 

z  (x,  0 )  =  zn(x) ,  %  (x,  0 )  =  yn(x ) 


(1.7) 


(1.8) 


z(x,T)=  z.( x ),  |r- (x,T)  *  y.  (x) 


(1.9) 

we  asked  if  there 


of  finite  energy,  i.  e.  C(Zq,  vQ  )  <  *  ,  C(Zj,  Vj)  <  •  ,  we  asked  if  there 
exists  ue  L2[0,T]  for  which  the  solution  of  (1.4),  (1.5),  (1.6)  correspond¬ 
ing  to  the  initial  state  (1.8)  assumes  the  desired  terminal  state  (1.9)  at  time 
t  =  T  .  The  answer,  a  qualified  "yes",  came  from  two  different  approaches  to 
the  problem.  The  relationship  between  these  two  approaches  has,  over  the 
years,  grown  ever  more  fundamental  and  has  led  to  a  great  many  very  interesting 
developments.  See  [34]  and  [45]  in  particular. 

The  first  method  explored  was,  as  we  have  already  indicated,  the  method 
of  characteristics.  If  we  let 

JW  (U0> 

and  consider  families  X+,  X”  of  "characteristic"  curves  satisfying 


.  v  ^ 

‘.\v  At  -^.1  A. 


+  C(x)  =  0.  -sg-  -  c(x)  =  0, 


(1.  11) 


respectively,  and  then  set 


v+(x,t)  =-|f-  (x, t )  +c(x)  ||  (x,t ), 
v‘(X,t)  =  ||  (x.t)  -  c(x)  ~  (x.t)  , 

we  see  readily  that  on  X+=  {(x+(t),  t)>  ,  X"  =  {(x-  (t),  t)},  respectively 
we  have 


v+(x  +  (t),  t)  =  c'(x+(t))  — 


(1.12) 


v~(x-  (t),  t)  =  c'(x-  (t))  — 


(1.13) 


t  =  t>  2Ti 


x  —  0 


t  =  T  <  2T. 


X'(L,T) 


X+(L,  0 


(x,t)  N 
\ 

\ 

*0  N 


x  =  0 


t  =  T  =  2T, 


Aj  x-(i;t) 


V(x,t) 

.  X+(L,0) 


Fig.  1.1:  The  Method  of  Characteristics 


Because  these  differential  equations  are  satisfied  on  different  families  of 
characteristics,  the  coupling  between  them  is  more  complicated  than  for  the 
usual  system  of  ordinary  differential  equations.  Nevertheless  there  is  a  method 
of  successive  approximations,  described  in  [30],  [36],  which  enables  solution 
of  these  equations  in  certain  regions  provided  with  appropriate  boundary  data. 
Such  a  region  is  the  roughly  triangular  domain  Aq  shown  in  Fig.  1. 1,  bounded 
by  t  =  0,  x  =  0  and  the  characteristic  X+(L,  0),  of  the  first  family 
described  by  (1.11),  passing  through  the  point  (L,  0) .  Together  with  the 
boundary  data  provided  by  (1. 5)  and  (1.  8),  it  may  be  seen  that  the  differential 
equations  (1.12),  (1.13)  determine  v+  and  v“,  and  hence  z(x,t),  through¬ 

out  the  domain  Aq  .  Similarly,  these  equations  together  with  the  data 
provided  by  (1.5)  and  (1.9)  determine  z(x,  t)  in  the  domain  Aj  bounded 
by  x  =  0,  t  =  T  and  the  characteristic  curve  X~(L,T),  described  by  the 
second  equation  in  (1.11)  and  passing  through  the  point  (L,T)  .  Thus  the 
initial  and  terminal  states,  described  by  (1.  8),  (1. 9),  together  with  the 
boundary  condition  (1.5)  determine  z(x,  t)  in  both  Aq  and  Aj . 

Whether  Afl  and  a^  are  disjoint,  or  have  a  region,  of  over¬ 

lap,  depends  on  the  time  T  allotted  for  control.  The  time  required  for  the 
curve  X+(0,  L)  to  pass  from  X  =  L  to  x  *  0  is 


and  this  is  also  the  time  required  for  X"(L,T)  to  pass  from  x  =  0  to  x  =  L. 
We  summarize  the  control  situation,  depending  on  the  relationship  between  T 
and  Tj  . 

Case  T  <  2T .  Here  AQ  and  A^  overlap  and  the  determinations  of 
z(x,t)  in  the  overlap  region  nQ=  aq  U  A^  provided  by  (1.8)  and  (1.9)  need 
not  and,  in  general,  will  not  agree.  There  can,  in  such  cases  of  disagreement, 
be  no  solution  of  (1.4),  equivalently  (1.12),  (1.13),  in  the  region 

Rj  =  {(x,  t)j0^x:£L,  OstsT).  The  control  function  u  (t ),  shown  in 
(1.  6),’  never  enters  the  picture  because  it  cannot  affect  the  solution  of  (1.  4)  in 
Aq  or  Aj  if  (1.8),  (1.9)  are  satisfied  at  t  =  0,  t  =  T,  respectively. 

Case  T  =  2T  .  Here  the  two  ’’domains  of  determ  inacy*',  Aq  and  a^, 
Just  fail  to  overlap;  their  boundaries  have  exactly  one  point  in  common,  t  =  T^, 
x  =  0.  The  Initial  and  terminal  conditions  (1.8)  and  (1.9)  determine  z(x,  t) 
in  Aq  and  A|  ,  respectively.  Another  process  of  integration  of  the  coupled 


*  .*«  V*  ,'*V‘  '• 

'  *  •  1  v*  v*  V*  V’  V  O 


differential  equations  (1.12)  and  (1.13)  permits  unique  extension  of  z(x,  t), 
equivalently  v+(x,  t),  v“(x,  t),  into  the  domain  J2  .  The  control  steering 
(1.  8)  to  (1. 9)  is  then  uniquely  determined  from  this  extension  and  (1.  6). 

The  determinations  of  z(x,t)  in  and  may  fail  to  match 
smoothly  at  the  point  p  :  x  =  0,  t  =  T.  .  This  results  in  discontinuities  of 

4-  «  * 

v  along  X  ( L,  0 )  and  of  v  along  X  ( L,  T  )  in  general. 

Case  T  >  2T .  The  only  difference  between  this  case  and  the  case 
T  =  2Tj  lies  in  the  line  segment  1:  x=0,  Tj  <  t  <  T  -  T^,  which  replaces 
the  point  p  of  the  case  T  =  2T^  .  Extension  of  z(x,t)  from  U  l 

into  n  cannot  be  carried  out  until  the  boundary  condition  (1.5),  which 
yields  9z/3t  ( 0,  t )  =  0  ,  is  augmented  by  arbitrary  data 

H  (0,t)  =  '  (°.t)  e  t  .  (1.  ' 

Once  this  is  done,  extension  of  z(x,  t)  into  0  proceeds  much  as  before. 
(See  [30],  [36]  for  details  of  the  extension  process. )  The  arbitrary 
function  £  (t )  can  be  designed  so  as  to  eliminate  discontinuities  of  the 
solution  along  X+(L,  0)  and  X“(L,T),  to  satisfy  some  criterion  of 
optimality  ( see  [  30]  e.  g.  )  or  to  fulfill  any  other  appropriate  design  objective. 

If  the  partial  differential  equation  (1.  4)  is  combined  with  boundary 
conditions  different  from  (1.  5),  (1.  6),  but  still  admissible  for  (1.  4),  the 

cases  T  <  2T^,  T  >  2T1  remain  as  above.  The  rather  delicate  situation 
at  T  =  2Tj  depends  on  the  specific  form  of  the  boundary  conditions.  For 
example,  the  boundary  conditions 

z  (0,  t )  =  0,  z  ( L,  t )  =  u(t) 

lead,  in  case  T  =  2T^,  to  a  situation  where  the  desired  control  is  not  unique; 
it  has  the  form 

u  ( t )  =  u(t)  +  yu(t) 

where  u(t)  is  a  non- zero  control  steering  the  zero  initial  state  into  the  zero 
final  state  and  y  Is  an  arbitrary  constant.  By  contrast,  the  boundary 
conditions 

f§(0.t)  =  0,  ||  (1,  t)  «  u(t)  (1.16) 

lead,  in  case  T  =  2Tp  to  a  situation  where  the  desired  control  u(t)  does 
not,  in  general,  exist.  (See  [31],  [37]  for  more  details. ) 

The  analysis  of  more  complicated  systems  of  hyperbolic  equations,  such 


as  the  Timoshenko  system  (1.  2),  (1.  3),  is  in  general  rather  complicated  but 
there  are  some  special  cases,  including  appropriate  boundary  conditions,  for 
which  the  analysis  is  fairly  simple.  In  [M]  a  discussion  is  given  permitting 
analysis  of  the  free  boundary  case 

||(0,t)  =  0,  +(0,t)-|£(0.t)=0,  (1.17) 

||  (L.t)  =  ux(t)  ,  4<(L,t)  -f£  (L,t)  =  u2(t)  .  (1.18) 

It  may  be  shown  that  all  cross -coupling  is  of  low  order  and  the  problem  is 
essentially  equivalent  to  two  problems  (1. 4)  with  boundary  conditions  (1. 16). 
Two  critical  times  are  involved.  With  (cf.  (1.14)) 

c  ,v)  -  ,«<*>  i  _  =  rl  dx 
cl(x)  ■  *  I p(x)  1  '  X1  J0  Cj(x) 

L 


(1. 19) 


*2<x>=  Tz  */ 


dx 


0  C2(X) 


(1.  20) 


'2'"'-'7M'  '  *2 

it  may  be  shown  that  finite  energy  states  are  controllable  if  and  only  if 

T  >:  2  max  {Tp  T2}  . 

The  essential  details  of  the  analysis  are  given  in  [30]  and  are  quite  similar 
to  what  we  have  briefly  outlined  here  for  (1. 4),  (1.  5),  (1.  6)  . 

It  is  immediately  clear  that  the  method  of  characteristics  is  specially 
adapted  to  controls  u(t)  acting  at  a  point,  as  in  (1.6).  This  is  true 
because  the  control  determination  occurs  at  the  very  last  stage  of  the  analysis, 
after  the  controlled  solution  has  been  computed.  If  the  control  u  (t ),  itself 
scalar,  acts  on  the  system  through  a  '’control  distribution  function"  g  (x),  as 
In  (cf.  (1.4)) 

a  9. 

(1.  21) 


'<*>  ^|-|j<S<X>ll>  =  9<*)u(t). 

ot 

homogeneous  boundary  conditions  (cf.  (1.5),  (1.6)) 

dz 


z(0.t)  =  0, 


dx 


(L,t)  =  0 


(1.  22) 


applying  at  the  boundaries,  we  face  what  appears  at  first  glance  to  be  a  rather 
different  situation  than  what  obtains  in  (1.  4),  (1.  5),  (1.  6),  for  even  the 


equations  corresponding  to  (1.12),  (1.13)  will  involve  the  unknown  control  u(t) 
in  this  situation;  one  cannot  proceed  by  filling  out  z  (x,  t)  in  successive 
domains  as  before;  a  completely  different  approach  is  required.  Such  an 
approach  can  be  found  in  the  study  of  moment  problems  -  a  technique  developed 


by  several  authors  (see  [3]  ,  [12]  ,  [15]  ,  [14]  ,  [23]  ).  The  technique  has 
the  advantage,  from  the  point  of  view  of  approximation  of  being  intimately 
connected  with  the  modal  representation  of  the  system  based  on  the  natural 
modes  of  vibration,  or  eigenfunctions  of  the  operator  -T{x)~^(d/bx)(s(x)(dz/bx)). 
It  is  known  (see  [1]  ,  [7]  )  that  the  operator 

Lz  =  -r(x)_1  (s(x)||  )  (1.23) 

with  boundary  conditions  conformable  with  (1.  22)  has  eigenvalues 

2  2 

Nc  =  (^TL  ]  “V“  +  ek  '  k  »  1.  2.  3,  .  • .  (1.  24) 

T1 

where  the  are  uniformly  bounded  and  is  related  to  c(x)  by  (1. 14). 

The  corresponding  eigenfunctions,  ^(x)  *  k  «  1,  2,  3,  •  •  • ,  form  an  ortho- 

normal  basis  for  L2  [0,  L]  (which  consists  of  the  same  functions  as 
2 

L  [0,  L]  but  has  the  inner  product 

(*»♦)*  J  r  (x)  <p  (x)  vJj  (x)  dx )  .  (1.25) 

0 

Every  finite  energy  solution  z(x,  t)  of  (1.21),  (1.22),  i.  e.  every  solution 
for  which  the  integral  (1.  7)  is  bounded  for  all  t,  can  be  expanded  In  the  form 


z(x,t)  =  Yj  zk(t>  ?k(x) 
k=l 

where,  if  we  assume  the  control  distribution  function  f(x)  has  the  expression 


g 


(x)  =  E  gk9k(x) 


k=l 


convergent  in  L  [0,  L],  the  z^(t)  satisfy 


+  Nc  zk  =  gk  u*tJ’  k  =  1,  2,  ••• 


(1.  26) 


Letting 


“k  =  N?  *  k  =  1,  2,  3,  •  •  • 


and  using  the  transformation 


1/iuv  -  l/l w. 


(1.  27) 


one  arrives  at  the  system 


It  may  be  seen  that  finite  energy  states  are  those  for  which 

Z  + VV2]  <  * 

k  =  1 

and  this  becomes,  in  terms  of  £k  » 

Z  (|iki2  +  kici2)  <  ■  • 

k=  1 

Integrating  (L  28),  we  have,  for  T  >  0  , 

iut-T  gv  -T  W(T-t) 

nvff)-e  K  nt(0)  =  -f  J  e  K  u(t)dt 


(1.  29) 


nkff)-e  nk(0)=  -f-  Jq  e 

iu,vT  g.  T  -iUkCT-t) 

CfcfD-e  Cfc( 0)  =  e 


u  ( t )  dt 


u  ( t )  dt 


Assuming  the  controllability  condition 

gfc  /  0,  k  *  1,  2, 3,  •  •• 

we  see  that  the  problem  of  steering  between  the  given  states  at  times  0  and  T 
reduces  to  the  moment  problem 


ZTiwTs 

f  e  *  f(s)ds  = 
J0 

.T  -iw-s 

/  e  K  f(s)ds  * 


\  * 


k  =  1,  2,  3, 


k  «  1.  2,  3, 


(1.  30) 


(1.  31) 


where  s  =  T-t ,  f  (s)  =  u  (T  -  s),  and 

i,  *  T  T 

ak  =  2(nk(T)  "  e  ^  V0,)*  Pk  =  2Uk(T)  -  e  *  Cic(°»  f1*32) 

are  square  summable. 

To  solve  the  moment  problem  we  resort  to  the  theory  of  nonharmonic 
Fourier  series  as  developed  by  Paley  and  Wiener  [24]  ,  Levinson  [  21]  , 
Schwartz  [42]  and  many  others.  (An  excellent  expository  treatment  [49]  by 
R.  Young  has  recently  appeared. )  The  following  is  known;  the  three  cases 
being  divided  in  a  manner  conformable  with  the  three  cases  discussed  earlier. 


iut.  t  -i  wtt 


Case  T  <  2T^ .  The  functions  e  ,  e 


k  =  1,2,3, 


linearly  dependent  in  L  [o,T]  in  a  rather  strong  sense.  Any  one  of  these 
functions,  indeed,  any  finite  number  of  them,  lie  in  the  closed  span  of  the 
remaining  functions  (which,  in  fact,  is  equal  to  the  whole  space  L  [0,  T]  ). 
As  a  result  the  moment  problem  (1.30),  (1.31)  cannot,  in  general,  be  solved. 


Suppose,  e.g. ,  all  but  finitely  many  of  the  Pk,  say  k=K+l,  K  +  2,  •••, 

were  equal  to  zero,  while  some  of  the  ,  3^,  k  =  1,  2,  •  •  • ,  K,  are  non-zero 
The  linear  dependence  just  referred  to  shows  that  such  a  problem  can  have  no 
solution;  the  equations 


T  iuvs 


f  e  K  f(s)ds  =  f  e  f(s)ds  =  0,  k  >  K 
0  0 

imply  that  the  same  equations  must  hold  for  k  s  K  . 


i<j  ^s 


-iu  vs 

e  *  .  k  =  1,2,3, 


Case  T  =  Tj.  Here  the  functions  e  ,  e  ,  k  =  1, 2, 3,  •  • 

form  a  Riesz  basis  for  L2[0,  ZT^]  .  Every  function  h  e  [0,  2T^]  has  the 
unique  convergent  expansion 


"  iuvs  -iu.  s 

h(s)  =  ]]  [hke  *  +  h_ke  K  ] 


(1.33) 


and  there  are  positive  numbers  c,  C,  such  that 

QO 

<='21|1>0Z.  *  £<|hk|2  +  |h.lt|2=SC2||h||2  (1.34) 

L  [0,2T]  *  *  L2[0, 2T] 

Further,  there  is  a  unique  dual  basis  of  biorthogonal  elements  p^,  P  k  e 
L2[  0,  2Tj]  such  that 

^1  ks  k  =  ±1,  ±2,  ••• 

J  e  P.(s)ds  *  5.  .  , 

0  1  K’1  l  -  ±1,  ±2,  . . . 


r  i.  k  =  t 

l  0,  k  /  i 


(1.  35) 


which  engenders  expansions  similar  to  (1.33),  the  roles  of  c,  C  in  the 
inequalities  parallelling  (1.  34)  being  reversed.  The  formal  solution  of  (1.  30), 
(1.  31)  is  then  uniquely  given  by 


f(s)  =  Y  [avPk(s)  +  \P_k  * 


(1.36) 


If  we  have 


Urn  |  g.  |  =  0  , 
k  —  ®  * 


(1.37) 


as  would  be  the  case,  e.g. ,  if  g  e  IT[  0,  L]  ,  then  the  conditions  for 
convergence  of  (1.  36)  are  more  stringent  than  just  the  square  summability  of 
the  given  by  (1.  32).  We  need 


As  a  consequence  we  can  steer  (1. 21)  from  any  finite  energy  initial  state  to  a 
dense  (in  the  energy  norm)  sub  space  of  final  states,  or  vice  versa,  but  we 
cannot  steer  between  arbitrary  finite  energy  states  during  [0,  2T^]  if  (1.  37)  is 
true.  The  case  of  boundary  control  (1.  6),  already  treated  by  the  method  of 
characteristics,  will  be  discussed  more  extensively  below.  In  that  case  the 
coefficients  gfc  In  (1.  26)  are  bounded  and  bounded  below.  The  result  is, 
in  that  situation,  that  we  obtain  the  same  result  this  way  as  by  the  method  of 
characteristics  -  given  finite  energy  initial  and  terminal  states,  there  is  a 
unique  control  u  e  L  [0,  21.  ]  steering  the  one  to  the  other. 


Case  T  >  2Tj  . 


The  main  difference  between  this  case  and  the  preced- 
-iu»i 


lwks 


**Tcs 

2. 


k  =  1,  2,  3, 


form  a 


ing  is  that  here  the  functions  e 

Riesz  basis  for  a  proper  subspace,  E,  of  L^fO,!]  .  The  biorthogonal 
functions  p^ ,  p  ^  exist,  but  are  unique  only  if  we  require  that  they  lie  in 
E  —  or  we  impose  some  comparable  condition.  If  we  agree  that  p^,  p  ^ 
belong  to  E ,  then  any  elements 


pk=pk+V  p-k  *  p-k  +  q-k 

with  qfc,  q_k  e  E1  C  L2[0,T]  still  form  a  biorthogonal  set  relative  to  the 

WS  -iubS  N 

e  ,  e  .  The  convergence  properties  of  series  involving  the  p^ , 

p_k  ate  much  the  same  as  in  the  preceding  case.  As  a  result  we  have  the 

same  control  capability  as  in  the  case  T  =  2T^  but  controls  are  not  unique  . 

Indeed,  if  u  is  a  control  steering  between  two  given  states,  the  family  of 

controls  u  +  r»,  ix  €  E*,  all  realize  the  same  control  objective.  Again, 

this  non- uniqueness  should  be  compared  with  the  similar  property  observed  for 

T  >  2Tj  in  applying  the  method  of  characteristics. 

Using  the  theory  of  distributions  and  related  material,  boundary  value 
control  situations  such  as  (1.  6)  can  be  included  in  the  same  framework  as 
(1.  26)  but  with  g  in  a  larger  space  than  L  [0,  L]  ;  g  should  be  a  linear 
functional  (in  general  unbounded  on  L  [0,  L]  )  whose  domain  includes  the 
domain,  jft(L),  of  the  self  adjoint  operator  L,  given  by  (1.23),  with  the 
given  homogeneous  boundary  conditions.  The  g^  are  the  values  which  g 
assumes  at  the  eigenfunctions  <p^  e  6  (L)  .  A  detailed  study  of  these 
"admissible  input  elements"  is  provided  in  [17]  .  In  this  way  a  unification 
of  the  boundary  and  distributed  control  cases  may  be  achieved.  One  consequence 
of  this  is  drat  the  biorthogonal  functions  p^,  p_^  which  play  such  an 


important  role  in  the  method  based  on  the  moment  problem  (1.  30),  (1.  31)  can 

actually  be  obtained  through  the  more  constructive  method  of  characteristics  as 
controls  steering  from  a  zero  initial  state  (say)  to  final  states  constructed 
using  a  single  eigenfunction  of  L  . 

We  began  our  discussion  here  with  the  Euler  beam  equation  (1. 1)  .  For 
definiteness,  let  us  add  a  distributed  control  term  (scalar  input)  and  specific 
boundary  conditions  so  that  we  have 


p(x) 

a2w 

+  -^T  (EI  (x) 

)  s  g(x)u(t) 

(1.  38) 

at2 

8x 

ax 

92w 

(0,t) 

_  a3w 

(0,t)  = 

0, 

(1.  39) 

ax2 

9x3 

2 

(L.t) 

a3w 

(L,t)  = 

0  . 

(1.  40) 

ax£ 

In  1969  I.  P.  Quinn,  in  his  doctoral  thesis  [27],  studied  the  controllability 
properties  of  a  class  of  systems  including  this  one.  Here  the  operator 


Aw  = 


p(x) 


(EI(x)-^) 


jT 

ax“  8x‘ 

on  the  domain  in  H4[  0,  L]  consisting  of  functions  obeying  boundary 
conditions  conformable  with  (1.39),  (1.40)  has  eigenfunctions  <p^(x)  forming 
an  orthonormal  basis  for  L2  [  0,  L  ]  and  the  corresponding  eigenvalues  X. 

4  p  1.  K 

grow  like  k  as  k  —  •  .  With  w  k  =  we  obtaln  a  system  similar 

to  (128),  using  a  transformation  like  (1.  27)  applied  to  the  second  order 

differential  equations  resulting  from  the  eigenfunction  decomposition: 


+\wk  =  gku(t),  k  e  0,  1,  2,  3,  ... 

( a  slight  modification  of  (1.27),  (1.28)  is  necessary  for  \Q  -  0  ;  see  [34]). 
Again  there  results 


and  the  energy  expression  for  (1.  38)  is,  equivalently, 
*/L[p(x)  (|^)2+EI(x)  (l^)2]dx 


u(t) 


(1.  41) 


or 


*  Z  t<«k>2+  W2i 

k=0 


j  £  inkr  +  kkr. 

k=0 

all  <  •  for  "finite  energy1 '  states.  4  _  ,  _ 

iw^S  -Iw^S 

Quinn  was  able  to  show  in  this  case  that  the  functions  e  ,  e  , 

lUlcS  •ig]/S 

k  =  0, 1, 2,  3,  •••  (for  k  =  0  replace  e  ,  e  by  1,  s)  are 

2 

linearly  independent  in  L  [0,T]  for  every  T  >  0  (this  result  by  itself 
had  already  been  obtained  much  earlier  by  Ingham  [18]  who  shows ,  in  effect , 
that  these  functions  form  a  Rlesz  basis  for  a  closed  subspace  of  L  [0,T]  for 
every  T  >  0 )  and,  additionally,  that  there  is  a  positive  number,  M(T  ), 
such  that  if  the  (non -unique)  biorthogonal  functions  p^fs),  p  Js)  are 
appropriately  selected  in  LZ[0,T]  ,  these  functions  are  continuous  and 
satisfy  the  pointwise  bounds 


|pk(s)|  *  M(T),  |p-k(s)|  =S  M(T),  s  e  [0.T]  . 


The  fact  that  the 


‘“ks 


(1.42) 


“i(UkS 


of  l/[0.T]  ,  T  >  0, 


e  form  a  Riesz  basis  for  a  closed  subspace 

Implies  that  initial  states  and  terminal  states  with 
(in  terms  of  (1.41))  expansion  coefficients  Hk  Q  ,  £k  Q  311(1  Tlk  j,  £k  j 
can  be  steered,  one  to  the  other,  during  [0,T]  with  ue  L2[0,t], 
provided  that 

•  2  Cv  n  2  *  \  ,  2  Ck  l  2 

En-5^-1  ><- •  ZU-5H-1  +  i-£ii><-. 

k=0  k  k  k=  0  k  k 

The  boundedness  property  (1.42)  shows  we  can  also  control  states  for  which 


Z  «■ 


k=0 


0» 

Z  (I 

k= 0 


I)  <  «, 


this  being  possible  with  a  control  function  u  (t)  uniformly  bounded  ana 
continuous  on  [0,T]  . 

We  have  noted  in  connection  with  the  Timoshenko  beam  system  (1.  2),  (1. 13) 


(1. 17),  (1. 18),  that  an  adequate  control  theory,  based  on  the  method  of  character¬ 
istics,  exists  when  we  have  two  separate  control  functions,  u^(t)  and  u2(t), 
with  which  to  control  the  lateral  deflection  and  shear  deformation  separately. 

An  open  question  is  the  adequacy  of  control  using  a  single  control  input,  so 
that  (1.17),  (1.18)  becomes,  e.  g., 

J^(L,t)  =  au(t),  *(M)-  j£  (L,t)  =  0u(t) 

z  z 

with  a  +  P  0  .  This  problem  is  a  special  case  of  the  more  general 
question  of  the  controllability  of  linear  hyperbolic  systems  of  dimension  n  =  2m, 
Involving  m  pairs  of  characteristics,  each  pair  describing  a  given  wave  mode 
propagating  in  two  opposite  directions,  by  means  of  fewer  than  m  control 
inputs.  Some  work  has  been  done  in  this  direction  by  R.  G.  Teglas  in  his 
thesis  [45]  and  by  N.  Wick  [47]  ,  but  it  is  safe  to  say  that  no  very  general 
criteria  for  this  problem  have  yet  appeared.  Particularly  valuable,  it  seems 
to  this  author,  would  be  a  study  of  the  Timoshenko  beam  system  from  the 
singular  perturbation  standpoint,  elucidating  the  behavior  of  solutions  and 
controllability  properties  as  the  modulus  of  elasticity  in  shear,  k  (x)  in  (1.  2), 
(1.  3),  tends  to  infinity. 

STABILIZATION,  CANONICAL  FORMS,  EIGENVALUE 
PLACE  ME  NT,  etc. 

As  all  practicing  engineers  will  know,  controllability  in  itself  is  rarely 
the  prime  goal  of  control  system  design.  Stability,  and  related  criteria  such  as 
robustness,  insensitivity  to  particular  input  frequency  bands,  etc. ,  are  more 
commonly  uppermost  in  mind.  Additionally,  there  is  the  question  of  state 
estimation  from  lower  dimensional,  noisy  observations  in  order  to  implement 
linear  feedback  control  policies.  These  subjects  have  been  pursued  almost  ad 
nauseum  for  linear,  finite  dimensional  systems.  In  the  case  of  distributed 
parameter  systems,  and  hyperbolic  systems  in  particular,  the  literature  on  this 
subject  remains  rather  sparse  and  spotty  in  its  coverage. 

As  in  the  case  of  linear  finite  dimensional  systems,  stability  and  stabili¬ 
zation  studies  for  linear  partial  differential  equations  have  tended  to  cluster 
around  two  dominant  approaches:  the  Liapounov  approach,  primarily  carried  out 
in  connection  with  systems  involving  some  form  of  "conservation  of  energy"  law, 
and  the  spectral  approach,  determining  if,  or  making  certain  that,  the 


eigenvalues  of  the  system  lie  in  an  appropriate  subset  of  the  left  half  plane.  The 
spectral  approach  suffers  from  the  disadvantages  of  greater  intricacy  of 
computation  and  the  need  to  show  that  the  spectrum  location  does,  in  fact, 
determine  the  asymptotic  behavior  of  the  system.  The  latter  brings  in  questions 
of  completeness  and  linear  independence  of  the  eigenvectors  of  the  system. 

We  will  begin  with  a  short  discussion  of  what  has  been  done  with 
Liapounov  methods.  On  the  theoretical  side  one  can  start  with  a  system 

X  S  Cx,  (2.1) 

C  generating  a  strongly  continuous  semigr  \p  S(t)  in  the  Hilbert  space  X 
(we  may  have  started  with  a  control  system  x  =  Ax  +  Bu,  set  u  =  Kx,  then 
C  =  A  +  BK )  .  We  set  up  a  quadratic  functional 

V (x)  =  (x,  Qx)  , 

where  Q  is  a  bounded,  positive,  self  adjoint  operator  on  X  with  Q  a  ql 
for  some  q  >  0 ,  to  serve  as  a  liapounov  function.  One  may  then  show  that 

for  t2  >  t|  and  x(t)  =  S  (t)xQ  a  "solution"  of  (2.1),  that 

(x(t2),  Qx(t2))  -  (x(tj),  Qx(tj)) 

^2 

*-/  (x(s),  Wx(s))ds 

*i 

for  some  positive  self  adjoint  operator  W  so  that,  in  some  sense  which  one 
needs  to  make  precise  in  individual  cases, 

C*Q  +  QC  +  W  =  0  .  (2.2) 

the  Liapounov  operator  equation,  is  satisfied.  An  important  result,  due  to 
Datko  [10],  states  that  if 

f"(x(t).  Qx  (t))dt  <  -  (2.3) 

0 

for  every  initial  state  xQ  e  X,  then  the  semigroup  S(t)  is  exponentially 
damped,  i.  e. 

US(t)U  s  Me"Yt  ,  t  a  0  , 

for  positive  numbers  M,  y  .  The  condition  (2.3)  is  satisfied  if  W>  wl 
for  some  w  >  0  ,  as  may  easily  be  verified. 

2 

Consider  the  linear  symmetric  hyperbolic  system  in  L*[0,  = 

(tZ[0,L])n 

E(x)|jr  +B(x)w  +t(x,t)  (2.4) 

where  E  (x),  A(x),  B  (x)  are  continuously  differentiable  mXm  matrices 
defined  for  x  e  [0,  L]  ,  E  (x)  symmetric  and  positive  definite,  A(x) 


symmetric.  The  wave  and  Timoshenko  equations  can  be  written  in  this  form. 

The  ’'energy"  usually  is  expressed  as 

i  rL 

e(t)  =  if  (w(x,t),  E  (x)w(x,t))dx  . 

0 

With  appropriately  "conservative"  or  "dissipative"  boundary  conditions  at 

x  =  0,  x  *s  L  ,  one  finds  that  for  t,  >  t. 

L  L  1  • 

e(t2)  -  CCtx)  =  /  (if  (w(x,t),  [B(x)  +  B(xr-A(x)]  w(x,t))dx 

tj  0 

,1* 

+  /  (w(x,t),  f (x,t))dx  >  dt  . 

0 

#  " 

If  B  (x)  +  B  (x)  -  A  (x)  is  uniformly  negative  definite  or  if  the  n  dimensional 
control  function  f  (x,  t )  may  be  arbitrarily  specified  as  a  function  of  x  and 
t,  one  may  use  feedback 

f  (x,t)  =  K(x)w(x,t)  (2.5) 

in  such  a  way  that 

*2  L 

C(t2)  -  ft(t,)  =  -/  /  (w(x,t),  W(x)  w(x,t))dx  dt  (2.6) 

^  ° 

with  W(x)  uniformly  positive  definite  and  symmetric  on  [0,  L]  .  Then  one 
can  apply  Datko1  s  result,  or  more  simple  arguments,  to  show  that  solutions  of 
(2.3),  (2.4)  are  uniformly  exponentially  damped  in  L^[0,  L]  norm. 

Note,  however,  that  if  B  (x)  +  B  (x)*  -  A  (x)  =  0  or  for  some  other 
reason  falls  to  be  positive  definite,  and  if 

f  (x,t)  =  D(x)u(x,  t)  or  f  (x,t)  =  D(x)u(t) 
with  dim  u(x,  t)  =  r<R  in  the  first  instance,  u  a  function  of  t  only  in 
the  second  instance,  then  we  cannot,  in  general,  achieve  (2.5)  with  W(x) 
uniformly  positive  definite.  Comparable  difficulties  arise  when  boundary 
control  is  employed.  In  such  cases  it  is  a  form  of  the  La  Salle  "invariance 
principle”  ( see,  e.  g.  [  19  ]  )  which  must  be  appealed  to,  rather  than  the  basic 
Liapounov  theory,  for  an  analysis  of  presumed  asymptotic  stability  properties  of 
the  system.  This  has  been  discussed  in  some  detail  in  [36]  and  [33]  and 
we  give  only  the  briefest  outline  here. 

The  "invariance  principle",  as  it  applies  to  finite  dimensional  systems, 
relies  heavily  on  the  compactness  of  the  "  u-limit  set"  of  the  system  in  order  to 
reach  the  final  conclusion  of  asymptotic  stability.  Comparable  compactness 
properties  associated  with  the  solutions  of  an  infinite  dimensional  system  are 


generally  difficult  to  realize  but  the  initial  attempts  to  extend  the  theory  nonethe¬ 
less  relied  on  establishing  some  sort  of  compactness  property.  One  of  the  first 
contributions  in  this  direction  was  due  to  Dafermos  [  9  ]  who  studied  weak 
damping  of  the  wave  equation,  relying  on  the  almost  periodic  nature  of  the 
system  solutions  to  provide  the  required  compactness.  Slemrod  [43]  studied 
the  boundary  damped  wave  equation  by  introducing  suitably  weakened  topologies 
—  as  compared  with  the  usual  topology  associated  with  the  energy  norm  —  and 
was  able  to  conclude  a  correspondingly  weakened  form  of  asymptotic  stability. 

Knowing  that  controllability  implies  stabilizability  in  the  case  of  autono¬ 
mous  finite  dimensional  linear  systems,  we  are  not  surprised  to  find  control¬ 
lability  playing  a  role  in  the  study  of  asymptotic  stability  and  stabilization 
properties  of  autonomous  infinite  dimensional  linear  systems.  This  is  discussed 
in  some  detail  in  the  paper  [28]  by  J.  P.  Quinn  and  the  author  and  also  in 
[33]  .  Systems  of  the  form  (2. 3),  but  with  the  control  appearing  in  the  bound¬ 
ary  conditions,  are  studied  in  [28]  prior  to  the  main  discussion  on  the  boundary 
damped,  higher  dimensional  wave  equation.  We  can  give  an  idea  of  the  flavor 
of  the  arguments  employed  using  a  simple  example  based  on  the  wave  equation 
(1.  4)  with  control  appearing  in  the  boundary  conditions  (1. 5),  (1.  6).  If  in  this 
system  one  employs  the  feedback  law 

u(t)  =  -y  (L,t),  (2.7) 

the  closed  loop  system  is  (1.4),  (1.5)  together  with  the  "closed-loop"  Robin 

type  boundary  condition 

(L,t)  +  y^(l,t)  =  0  .  (2.8) 

Here  a  short  computation  shows  that  with  the  energy  e(t )  defined  by  the 
expression  (1.7)  we  have,  for 


t2  >  tj. 


e(t2)  -  eftp  =  -YP  (I)/  ^(L,trdt.  (2.9) 

It  is  not  feasible  to  fit  this  situation  into  the  general  pattern  based  on  the 
Liapounov  operator  equation  (2.  2)  but,  since  we  expect  (correctly)  that,  along 
with  (1.4),  (1.5),  (2.8) 


^  ( L,  t )  2  0  =$>  w(x,  t)  2  0  , 

an  "invariance  principle"  type  of  argument  appears  to  be  in  order.  But  we 
will  use  a  variation  on  this  procedure  which  makes  use  of  the  controllability 
already  established  in  Section  1.  Let  v  (x,t)  be  a  controlled  solution  of 


(1. 4),  (1. 5),  (1.  6),  u  (t)  being  selected  so  as  to  steer  the  initial  state 

v  (x,  0)  =  w  (x,  0),  (x,  0)  =  jjjr  (x,0),  (2.10) 

agreeing  with  the  initial  state  of  the  solution  w  (x,  t)  of  (1.4),  (1.5),  (2.8), 
to  the  zero  final  state 

v(x,Tx)  =  0,  ^(x,^)  =  0.  (2.11) 

Tj  as  described  earlier.  Defining  the  "energy  inner  product" 

(w(.,  t),v(*  ,  t))  =  f  [p  (x)  (x.t)  |f  (x,t)  +  p(x)  |^(x,t)|j(x,t)]  dx 

it  is  found,  using  (1.  5),  (2. 8),  (2. 10),  (2. 11),  that 

<w(.,0),  v(  • ,  0))  -<w(..2T),  v(.,  2T,)>  =  jjw(*,  0)|J2  = 

2TL  1  1  c 

=  -P(L)  /  [&(L,t)|f(L,t)+-^(L,t)^-(L,t)]dt 


=  P  (L>/~  1  IF  tv|r  (L,t)  +  |f(L,t)]dt 


-P«J  1  ^(L.t)[y  |f  (L.t)  +u(t)]dt. 

0 

Here  ||w(-  ,  0)|J^,  the  energy  norm  at  t  =  0,  is  2e(0)  .  Applying  the 
Schwartz  Inequality 


4e(0)2  *  P  (L)  / 


1  3w 


(L,t)*dt 


(2. 12) 


•  ^  1  (Y^(I*ft)  +u(t))2dt. 

A  slightly  more  detailed  study  of  the  control  problem  for  (1. 4),  (1.  5),  (1.  6)  in 

the  case  T  =  2T.  (or  T  >  2T. )  shows  that  control  from  an  initial  state 
w(x,  0),  (x,  0)  to  0,0  at  time  2T^  is  realized  with  a  control  u(t) 

which  satisfies  2X 

/  u(t)2dt  s  K0  5(0) 

0  0 

and,  for  the  resulting  controlled  solution  we  have 

2T 

/  1  If  (L,t)2dt  ^  Kxe(0) 

for  certain  positive  constants  K_  and  K,  .  Then  (2. 12)  easily  yields 
2T  u  i 

YP fl*)/  ( L, t ) 2 dt  s  - -  =  Ke(0) 

0  dt  2  (Kq  +y  n.)  6(0) 


(2.13) 


and,  setting  t^  =  0  ,  t2  =  2T^  in  (2.9),  we  have 

e(2TL)  s  e(0)  -  Ke(o)  =  (l-K)  e(o)  . 

Since  e( 2Tj)  is,  from  (2.9),  (2.13),  positive  and  less  than  or  equal  to  £(0) 

we  conclude  0  <  1  -  K  <  1  . 

Repeating  the  above  argument  on  successive  intervals  [0,  2T^]  ,  [2Tj, 

4Tj]  •••  (2kTp  ZOc+lJTj]  , . . .  and  using  the  monotonicity  of  e(t),  as 
implied  by  (2.9),  we  conclude  that  e(t)  decays  exponentially  to  0  as  t  — «. 
The  same  general  argument  can  be  used  with  a  fairly  wide  class  of  boundary 
damped  linear  symmetric  hyperbolic  systems  (2.  4)  and  with  many  other  systems 
which  are  energy  conserving  in  the  uncontrolled  situation  and  suitably  strong 
controllability  properties.  The  Timoshenko  system  (1.  2),  (1. 3),  with  appropriate 
boundary  conditions,  is  in  this  class.  As  far  as  the  author  is  aware,  the  Euler 
beam  model  (1. 1)  has  not  yet  been  studied  from  this  point  of  view. 

The  spectral  approach,  as  we  have  already  indicated,  involves  a  direct 
analysis  of  the  eigenvalues  and  eigenfunctions  or,  more  generally,  the  spectrum 

and  invariant  subspaces,  of  the  generating  operator  C  for  a  given  system 

•  • 

x  =  Cx,  possibly  derived  from  a  control  system  x  =  Ax  +  Bu  by  the  use  of 
linear  feedback  u  =  Kx  so  that  C  =  A  +  BK  .  A  fairly  common  case,  which 
can  be  treated  with  minimal  difficulty,  arises  when  all  but  finitely  many  of  the 
eigenvalues  of  C  have  negative  real  parts.  Under  genetically  valid  control- 
lability-type  conditions  it  is  then  possible  to  move  the  unstable  eigenvalues 
into  the  left  half  plane  while  either  keeping  the  stable  eigenvalues  fixed  or  else 
maintaining  a  certain  margin  of  stability.  Work  of  this  sort  has  been  carried  out 
by  Trlggiani  [46]  ,  Sakawa  [40]  ,  [41]  and  others. 

A  somewhat  more  challenging  task  arises  when  one  starts  with  a  system 
having  infinitely  many  eigenvalues  in  the  closed  right  half  plane  (usually  one 
considers  a  conservative  system  wherein  all  of  the  eigenvalues  of  C  are  purely 
imaginary)  and  one  attempts  to  devise  a  feedback  law  to  move  all  of  these 
eigenvalues  over  into  the  open  left  half  plane.  A  number  of  procedures  have  been 
examined  in  this  connection. 

In  [32]  a  second  order  system  with  scalar  control 

x  +  Ax  =  bu,  x,  b  e  X  ,  (2. 14) 

is  studied,  X  being  a  real  Hilbert  space  and  A  an  unbounded  positive  self 
adjoint  operator  on  X  .  Assuming  that  A  has  a  Riesz  basis  of  eigenvect  r's 
k=  1, 2, 3,«**,  in  X,  and  corresponding  positive  eigenvalues 


(2.15) 


increasing  with  k  ,  k  =  1, 2,  3,  •  •  • ,  x  and  b  may  be  expanded  as 

*=  i  *** k -  b  ■  2  bk4’ 

k=  1  k=l 

convergent  in  X  ,  with  square  summable  coefficients.  We  assume  the 
minimal  condition  for  approximate  contro liability 

\  /  0  *  k  =  1,2,3,  •••  . 

•  • 

The  energy  form  is  £  [  (x,x )  +  (x,  Ax )]  =  j  and  elementary  computations 
show  that  for  (2. 14)  and  for  any  T  >  0 

rT  * 

e(T  )  -  e(0)  =  /  (x(t),b)u(t)dt  .  (2.16) 

0 

It  follows  that  with 


u{t)  =  -y(x(t),b)  (2.17) 

the  energy  £(t)  is  non- increasing  with  increasing  t  .  So  far  this  is 
basically  a  Liapounov  approach  employing  what  is  known  in  the  engineering 
literature  as  an  ILAF  (Identical  Location  of  Accelerometer  and  Forces) 
approach.  The  resulting  closed  loop  system  is,  still  in  second  order  form, 

x  +  Bx  +  Ax  s  o  (2. 18) 

with  B  defined  by 

Bx  =  y(x,b)b.  (2.19) 

With  y  =  x ,  one  may  consider  the  equivalent  first  order  system  in  XX  X, 

(*)-  c{y)'c‘(*  -b)  ■  <2-2o: 


and  ask:  what  are  the  eigenvalues  and  eigenvectors  of  C?  It  is  here  that 
one  leaves  the  second  method  of  Liapounov  and  returns  to  his  first.  In  [32] 
a  perturbation  analysis  is  carried  out,  valid  for  small  values  of  y  in  (2.17), 
(2.19).  It  Is  shown  that,  under  the  separation  assumption 

wk+l  "  wk  “  d  >  0’  =  •  (2. 21) 


(2.  21) 


the  eigenvalues  of  C,  which  for  y  =  0  are  ±iwk,  k=  1, 2, 3,*»*,  all 
have  negative  real  parts  for  y  >  0  and,  moreover,  designating  the  perturbed 
eigenvalues  by  £k(y).  k  =  ±  1,  ±  2,  ±3,---,  £k(0)  =  iwfc,  £_k(°)  =  -iwk, 
we  have  (cf.  (2. 21)) 


Ck(Y)  =  iwk  -  |  Jbjl2  +  ®(y2. 


),  k  - 


(2.  22 


),  k  —  »  .  (2.22) 


£_k(Y )  =  -iwk  -  |-|bk|2  +  (S{v2 


1 


(wvi‘ 


It  is  also  possible  to  show  that  the  perturbed  eigenvectors  continue  to  form  a 
Riesz  basis  for  the  space  XXX.  From  this  it  follows  that  all  solutions  of 
(2. 18)  tend  strongly  to  zero  in  the  energy  norm,  though  not  at  a  uniform  expon¬ 
ential  rate. 

Following  Wonham' s  initial  results  [48]  on  the  finite  dimensional  case,j 
there  has  been  considerable  interest  displayed  in  the  question  of  spectral  ! 

determination  via  linear  feedback  for  distributed  parameter  systems.  In  terms  ! 
of  the  system  (2. 14),  equivalently,  ! 


u. 


(2.  23) 


with  initial  (u  =  0)  eigenvalues  ±  iwk,  k  =  1,  2,  3,  •  *  • ,  the  question  may  j 

be  phrased  as  follows:  we  suppose  use  of  a  linear  feedback  functional  \ 

1 

u  =  (A2x,  kx)  +  (y,k2),  kr  k2  e  X,  (2.  24)  j 

bounded  relative  to  the  energy  norm 
XXX.  With 

KjX  =:  (A^x,  kj)b, 
the  closed  loop  system  is 


(x,Ax)  +  (y,y)  =  (A2x,  A2x)  +  (y,  y)  in 
K2y=(y,k2)b  (2.25) 


(2.  26) 


One  can  now  ask:  What  eigenvalues  can  be  achieved  for  the  closed  loop  I 

system  (2.24)  by  appropriate  selection  of  k^  k2  in  (2.25)?  For  some  j 
time  the  author  was  under  the  impression  that  his  approach  via  canonical  forms  j 
[35]  (more  on  this  below)  was  the  first  treatment  of  this  question  but,  in  i 
fact,  it  appears  that  this  credit  must  go  to  Prof.  Sun  S.  -H.  of  Szechuan  I 

f 

University  who  treated  this  problem  by  a  more  sophisticated  application  of  the  ; 
perturbation  technique  used  by  the  author  in  [32]  to  obtain  the  result . 

Sun  was  able  to  show,  with  an  assumption  similar  to  (2.  21)  and  the  Riesz 
basis  assumption  on  the  open  loop  eigenvectors,  that  the  totality  of  spectra,  Yjt 
achievable  by  use  of  (2.25)  coincides  with  sequences  £k,  £  k,  k  = 

2,  3,  •••  for  which,  assuming  the  ^  ^  0  as  before. 


His  very  important  paper  has  been  translated  by  Ho  L.  -F.  in  [44]  .  Some 
comparable,  but  necessarily  weaker,  results  have  been  obtained  by  Reid  in  his 
thesis  [29]  for  the  equation  of  linear  surface  waves  where  (2.21)  is  not 
satisfied  and,  in  fact,  11  m  (wk+1  -  w^)  =  0  .  Other  results  in  this 

direction,  for  hyperbolic  systems  of  various  types,  have  been  obtained  by  Clark 
[5]  ,  [6]  and  by  Ho  in  his  thesis  [16]  . 

Much  of  the  initial  impetus  for  the  study  of  control  canonical  forms, 
both  for  finite  and  infinite  dimensional  systems,  came  from  the  spectral 
determination  question  discussed  above,  but  the  subject  is  interesting  in  its 
own  right  and  shows  some  promise  of  being  adaptable  for  "real  world"  control 
implementation.  The  reader  will  recall  that  a  finite  dimensional  controllable 
system 

x  =  Ax  +  bu,  x  e  Rn , 

with  scalar  control  u  is  equivalent,  via  a  state  space  similarity  transform¬ 
ation  (see  [20],  [35])  to  a  system  in  rational  canonical  form  corresponding 
to  the  n-th  order  scalar  equation 


y(n>  +  a.yf"-1’  + 


y  +v= u- 


(2.  27) 


where 


p(X)  =  det(XI-A)  =  kn  +  ax  \n_1  + 


a  +  a 
n-l  n 


is  the  characteristic  polynomial  of  the  matrix  A.  Comparable,  but  somewhat  I 
more  intricate,  results  are  available  for  systems  with  higher  control  dimension  ! 

t 

[20]  ,  [2]  .  In  [38]  we  note  that  if  one  employs  a  scalar  linear  observation  j 

I 

y  =  h*x  =  (x,  h),  (2.  28)  j 

there  is  exactly  one  observation  vector  h  e  Rn  for  which  (2.  28)  satisfies  ! 
(2.  27);  for  general  h  the  right  hand  side  will  involve  the  derivatives  of  u  1 
of  order  at  u  -  1  .  Systems  (2.  27)  are  particularly  easy  to  deal  with. 

Closed  loop  eigenvalues  ^ may  be  realized  simply  by  forming 
the  polynomial 


q(M  =  TT  (x-<k)  =  >^n  +c, 

k  =  1  K  1 


\n -1  + 


+  c  .X.  +  c 
n-l  n 


and  determining  u  by  linear  feedback  on  the  observation 


and  its 


derivatives , 


n 


(n-k) 


u  =  l  <ate-ck)y 

k  =  1  ! 

Apparently  less  well  known,  but  quite  obvious,  is  that  the  control  problem  for  i 

i 

(2.  27)  is,  in  a  sense,  trivial.  Let  us  suppose  the  initial  instant  is  taken  to  i 
be  t  =  0  and  control  is  to  be  effected  during  0  s  t  s  T  .  Let  the  initial! 
state  be  specified  by 
r(n-k) 


k  =  1,  2, 


k  =  1,  2, 


y . <°>  ■  yn-k+l  • 

and  the  terminal  state  by 

y(n-k) «.%  y 
y  n-k+1  ' 

If  y(t)  satisfies  (2.29)  and 

y*n)(t)  =  v(t),  OstsT 
then  we  see  readily  that  for  k  =  1,  2,  •  •  • ,  n 

.  .  '  k  v-f  t 

,(n-k)  v  t 


,  n 


,  n 


y-  '  (t)  =  l  yn_(+1  +  /  (t -S  )k_1  V(s)  ds 

£  =  1  U 


and  (2.  30)  is  achieved  just  in  case 

T  ,k-£  k 


,k-£ 


/  (t-s)K-,v(s)ds  =  Yn.k.1-2  yn.t_!  -^rrn 

V  *  —  1 

k  =  1,  2, 


1=1 


(2.  29) 


(2.30) 


(2.  31) 


,  n  . 


This  is  easily  solved  for  v  in  various  function  classes,  e.g.  polynomials  of 

1 

degree  s  n  -  1,  etc.  and,  it  should  be  noted,  the  solution  has  nothing  to  do  i 
with  the  coefficients  in  (2.  27)  so  the  calculation  can  be  carried  out  once  for 
any  given  T  and  recorded  for  use  ever  after.  Then  in  a  given  canonical 
system  (2.27)  we  need  only  set 

n  j 

u  (t)  =  v  (t)  -  £  a^^tt)  (2.  32)  i 

k=  1  | 

to  realize  the  desired  control  objective.  ’ 

Since,  in  a  given  control  context,  it  is  not  likely  that  the  available 
observation  (2.  28)  will  be  the  particular  one  for  which  (2.  27)  obtains,  the 
above  result  might  seems  to  be  a  generally  useless  curiosity.  It  turns  out, 
however,  that  in  canonical  form  theory  there  is  a  counterpart  to  the  more  widely 
known  observer  theory.  If  C  is  any  nXn  matrix  whose  minimal  and 


characteristic  polynomials  coincide,  it  is  possible  to  select  (non-uniquely)  r, 
d  and  j  such  that  the  augmented  system 
x  =  Ax  +  bu 


z  =  ry  +  Cz  +  du  ( =  rh^x  +  Cz  +  du  since  y  =  h  x ) 


(2.  33) 
(2.  34) 


with  augmented  observation 

w  =  y  +  j ' z  =  h’x  +  j'z 


(2.35) 


is  in  canonical  form,  so  that  for  some  coefficients  o^,  a 2,  • 


’  “2m 


w^2n^  +  +  •  *  •  +  a2m-l  W'  +a2mW=  u 


The  adjoined  system  (2.  34)  can  be  realized  electronically,  just  as  an  observer, 
system  is,  and  the  considerable  freedom  in  choice  of  C  ,  r,  d  and  j 
provides  much  design  flexibility.  In  some  cases  the  dimension  of  (2.  34)  can  ; 
be  reduced.  The  proof  that  (2.33),  (2.34),  (2.  35)  can  be  made  a  canonical  I 

system  appears  in  [38]  . 

A  parallel  control  canonical  form  theory  has  been  developed  for  certain 
hyperbolic  distributed  parameter  systems,  involving  neutral  functional  equations! 
in  place  of  the  n-th  order  scalar  equation  (2.  27).  The  theory  is  quite  complex!, 
especially  as  it  applies  to  partial  differential  equations  with  variable 
coefficients  (see  [35],  [16],  [38],  [39]  e.g. ).  To  give  an  idea  how  the 
theory  is  developed  we  will  consider  the  constant  coefficient  case  of  (1.  4) 


which,  without  loss  of  generality,  we  can  take  to  be 


32w 


32w 


=  0,  t>0, 


1, 


(2.  36) 

at~  ax'* 

w(°,t)  =  °,  ~  (1,  t )  =  u  (t )  .  (2.37) 

The  normalized  eigenfunctions  of  the  corresponding  homogeneous  system  are 


rf.  (x)  =  ^Tz~  sin 


2k-l 


Setting 


wk  = 


k 
2k -1 


TTX,  k  =  1,  2,  3, 


(2.  38) 


7T 


and  forming  the  expansions 


W(X,t)  =  Yj  W]C(t)  » 

k=  1 


(2.  39) 


I r(x-t)=£ 

k=l 

followed  by  the  transformation 


(2.  40) 


(2.41)  1 


we  have,  for  k  =  1,  2,  3,  •  •  • 

\  *  lwknk  +  u}=-  u(t>-  <*  ■  -lwkCk  +  Ljj=~  “ 

Consider  now  the  neutral  delay  equation 


(2.  42) 


y(t  +  2)  +  y(t)  =  u(t+2)  .  (2.43) 

The  characteristic  function  of  the  homogeneous  equation  is 

p  (X)  =  e2^  +1  =  2e^  cosh  X 

and  the  zeros  of  p(X)  are  precisely  the  eigenvalues  ±  iw^  appearing  in 
(2. 42)  .  The  transfer  function  for  (2.  43)  is 

T  fXl  -  -  JL  s*nh  V  +1  (2  4^ 

T0(X)  “  2k"  ,  "  2  cosh  fc.  +  2  (  * 

e  +i 


L  sinh  X  1^ 

2  cosh  X  2 


(2. 44) 


which  can  be  rewritten  as 


Vx>=  l  T TT-r 

k=l  +  “k 


,  1  _  2k -1 

+  2  '  “k  “  2 


(2.44) 


If  we  define  an  observation  y  (t)  on  (2.  42)  by 

OB 

y(t)  =  E  thknk(t) +gk^k(t)]  +  iu<t> 

k  =  1 

the  transfer  function  for  y  is,  formally. 


■  ft  /  nk-l  — 

y  V-1* 

-1_  K-  l“k 


V-1^'1  7T 

X  +  i  Ut 


which  may  be  seen  to  agree  with  (2.  44)  just  in  case 

Using  (2.  38),  (2.  39),  (2.41),  (2.42)  it  may  be  seen  that  this  choice  of  h^,  g^J 

corresponds  to 

y(t)  =  \  +  j*(t)  =  \  (-^(l,t)  +  ~-(l,t))  .  (2.45)  i 

This  observation  on  (2.  42),  and  no  other,  satisfies  the  scalar  equation  (2. 43) 
which  serves  as  the  control  canonical  form  for  (2.42).  The  details  of  the  above 


calculations  and  some  idea  of  the  form  of  a  general  theory  appear  in  [38]  and 
[39]  . 

If  the  canonical  observation  (2.  45)  were  actually  available,  so  that 

we  have  (2.43),  its  usefulness  is  quite  clear.  For,  with  the  causal  feedback 

law  2 

u(t  +  2)  «  (1- Y)y(t)  -  /  c  (s)y  (t  +s)ds 

0 

(2.  43)  transforms  to 

2 

y  (t +2)  +  yy(t )  +  f  c(s)y(t  +  s)ds  =  0 

0  Ck1 

and  it  is  known  from  [35]  ,  [44]  that  the  exponential  solution  e  , 

of  (2.  47)  can  be  made  such  that 


(2. 46) 


(2.47) 

C-k* 


£k  =  i^k  +  “  +  ek*  C_k  =  iojk  +  a  +  e-k  ' 
where  a  is  a  complex  number  (ordinarily  negative)  determined  by  y  and 
are  arbitrary  complex  numbers,  determined  by  cc  L  [0,2],  such 

OB 

.2  .  I  ,2, 


ek’  e-k 
that 


Z  <Kr  +  l«.kl  >  <  -  • 

k  =  i 

It  may  be  shown  that  these  are  the  eigenvalues  of  the  closed  loop  system  (2.36), 
(2.37),  (2.45),  (2.46). 

In  a  given  application,  however,  it  is  entirely  likely  that  the  particular 
"canonical"  observation  (2.46)  will  not  be  available.  Indeed,  in  the  example 
indicated,  since  this  observation  is  taken  at  the  same  point  where  control  is 
applied  and  might,  therefore,  be  subject  to  a  certain  amount  of  noise 

disturbance,  it  might  not  be  desirable  to  use  this  observation  in  practice.  To  | 

» 

illustrate  the  use  of  the  technique  of  canonical  augmentation  (or  "canonical  j 
compensation",  perhaps)  let  us  consider  the  same  system  (2.3  6),  (2.37), 
but  suppose  the  available  observation  is  j 

(2.48)  j 


y(t)  =  (o, t) 


It  is  not  hard  to  show  in  this  case  that  y  (t)  satisfies 

y(t  +  2)  +  y  (t )  =  u(t  +  l) 


(2.  49) 


rather  than  (2.  43).  This  "central"  control  canonical  form  is  not  as  usable  as 

the  "backward"  form  (2.43)  because,  unlike  (2.46), 

,2 

u(t  +  l)  =  (l-y)y(t)  -  J  c (s)  y(t  +  s)ds 

0 


is  not  a  causal  feedback  law  and  cannot  be  implemented.  But  now  couple  (2. 49) 
with 

z(t+2)  +  Pz(t)  =  au(t+2)  +bu(t+l) 


+  cy  (t  +  1)  +  dy  (t) 


(2.50) 


and  let 


w(t)  =  y (t)  +  z(t)  . 


One  ordinarily  will  take  |  p  J  <  1  so  that  the  homogeneous  part  of  (2. 50)  is 
asymptotically  stable,  thus  avoiding  the  growth  of  parasitic  solutions  in  the 
compensator.  Since 


[y(t+4) +y  (t+2)  -  u(t+3)]  +  p  [y  (t+2  )  +  y  (t)  -  u  (t+1)]  =0 


while 


[z(t+4)+pz(t  +  2)  -  au  (t  +  4)  -  bu  (t +3)  -  cu(t +1)  -  cy(t +3) -dy  (t +2)] 

+  [  z(t  +2)  +  pz  (t)  -  au  (t  +  2)  -  bu  (t+1)  -  cy  (t  +1)  -  dy  (t)  ]  =0 
we  find  that 

w(t+4)  +  (l+  p)w(t+2)  +  pw(t)  =  au  (t  +  4)  +  [1+b]  u(t+3) 

+  au (t  +2)  +[p  +  b]  u (t  +1)  +  c[y(t+3)  +  y(t+l)] 

+  d[y(t  +  2)  +y(t)]  *  (using  (2.49)) 
au (t  +  4)  +  [1+b]  u  (t  +3)  +  [a+c]  u(t  +  2)  +  [p+b+d]  u(t+l)  . 
Then  it  is  easy  to  see  that  with 


a  =  1,  1  +  b  =  a  +  c  =  p  +  b+  d  =  0, 


i.  e.  with 


a  =  1,  b  =  -1,  c  =  -1,  d  =  1  -  p  , 
we  arrive  at  the  "backward*'  canonical  form  satisfied  by  w(t) 

w(t  +  4)  +(l  +  p)  w(t+2)  +  pw(t)  =  u(t  +  4) 

for  which  causal  feedback  laws 

u  ( t + 4)  =  -YjW(t+3)  +  [1+p  -y2]  w  (t  +2)  -  Y3  w(t  +  l) 

4 

+  [P“Y4]  w(t)  -  J  c(s)w(t  +  s)ds 
may  be  implemented,  yielding  overall  closed  loop  systems 


(2.51) 


(2.52) 


w(t  +  4)  +Y]W(t+3)  +y^w(t  +  Z)  +YjW(t+l)  +Y4w(t) 

/>  ^ 

+  I  c(s)w(t+s)ds  =  0. 

0 

It  is  necessary  to  check  separately  that  the  system  (2.49),  (2.50),  (2.51)  is  • 

observable  in  any  given  case. 

The  exponential  solutions  of  (2.  52),  and  hence  the  eigenvalues  of 
(2.36),  (2.37),  (2.48),  (2.50),  (2.51)  may  be  determined  with  the  same  ! 

flexibility  as  already  noted  for  (2.  47).  This  is  discussed  in  some  detail  in  the 
thesis  of  R.  G.  Teglas  [45].  A  complete  theory  of  canonical  compensation 
for  hyperbolic  systems  remains  to  be  developed  but,  we  hope,  the  example 
given  here  gives  reason  to  believe  that  the  method  is  a  promising  one.  It  is 

i 

clear  that  there  are  some  connections  with  observer  theory  as  developed  in  j 
[22]  and  elsewhere;  these  connections  remain  to  be  worked  out.  < 
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Admissible  Input  Clements 
for  Systems  in  Hilbert  Space 
and  a  Carleson  Measure  Criterion^ 


by  E.F.  Ho+  and  D.  L.  Russell++ 


Abstract 

We  study  the  control  system 

x  s  Ax  +  bu,  x  €  X  ,  u  scalar, 

where  A  generates  a  semigroup  on  the  Hilbert  space  X  but,  in  general, 
the  control  input  element  b  /  X,  Many  boundary  value  control  systems, 
point  control  force  situations,  etc. ,  can  be  studied  in  this  context  .  We 
define  and  analyze  "admissible"  input  elements  b  and  develop  sufficient 
conditions  for  b  to  be  admissible  in  terms  of  the  Carleson  measure  theorem 
of  Hp-theory  . 
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1.  Introduction . 


One  commonly  studies  linear,  time  invariant  control 


systems  in  a  Banach  space  X  in  the  form 

x  =  Ax  +  Bu,  x  €  X  ,  u  e  U , 


(1. 


where  A  is  the  generator  of  a  strongly  continuous  semigroup  of  bounded 
operators  {S  (t)|  t  2  0}  on  X  and  B  is  a  bounded  operator  from  the 
control  space ,  U,  into  X.  If  u:[0,®)-»  U  is  locally  (Bochner) 
integrable,  generalized  (or  "mild")  solutions  of  (1.1)  corresponding  to 
an  initial  state 


x  (0)  =  xQ  e  X 


can  be  represented  by  the  "variation  of  parameters"  formula  (see,  e.g.  [3] 


[11]) 

x(t)  =  S(t)xQ  +  /  S(t-  s)Bu(s)ds 


(L 


and  a  number  of  properties  of  x(t)  thereby  deduced  . 

It  is  well  known,  however,  that  most  of  the  "interesting"  infinite 
dimensional  control  systems  do  not  arise  this  way  because  the  degree  of 
controllability  of  a  system  (L  1)  with  B  bounded  is  rather  restricted  if,  as 
is  usually  the  case,  U  is  finite  dimensional  or  for  some  other  reason  the 
operator  B  is  compact  .  Indeed,  most  of  the  mathematically  intriguing 
examples  arise  in  the  context  of  partial  differential  equations  with  boundary 
value  control  inputs,  control  forces  exerted  at  isolated  points,  etc. ,  and  in 
the  context  of  functional  equations  which  involve  values  of  the  control  of 
discrete  instants  ,  viz. ;  u  (t) ,  u  (t  -  T^), ...  ,  u  ( t  -  Tn)  .  In  each  of  these 
cases  the  formulation  (1. 1)  is  inadequate  and  one  must  consider  input 


operators  '  B  whose  range  is  not  restricted  to  the  space  X  . 

A  number  of  authors  have  addressed  the  problem  of  interpretation  of  (1. 1) 
for  operators  B  of  rather  general  type  .  We  particularly  cite  the  contributions 
of  Curtain  and  Pritchard  [3]  ,  Zabczyk  [22]  ,  Fattorini  [6]  ,  and  Washburn  [20] 
It  seems  fair  to  say  that,  as  brought  out  in  [  3  ]  ,  the  theory  is  more  extensive 
and  generally  applicable  in  the  case  of  systems  of  "diffusion  type",  ordinarily 
involving  holomorphic  semigroups,  than  in  systems  of  "wave"  or  hyperbolic 
character  . 

In  the  present  article  we  shall  restrict  our  attention  to  spaces  X  which 

are  separable  Hilbert  spaces  and  to  finite  dimensional  control  spaces  U  . 

Taking  U  to  be  Rm  ,  (1. 1)  becomes 

m 

x  =  Ax  +  Yj  bj  (1.  2) 

J=1 

where  b^  is  the  ontrol  input  element  associated  with  the  j-th  control 
component  u^  .  Since  every  solution  of  (1.2)  is  a  linear  combination  of 
solutions  of  x  =  Ax  and  the  individual  systems  x  =  Ax  +  b^  u^  , 
j  =  1,  2, . . .,  m  ,  we  may,  without  loss  of  generality,  confine  our  discussion 
to  systems 

x  =  Ax  +  bu  (1.  3) 

wherein  the  control  u  is  scalar  valued.  Much  of  our  theory  can  be 
extended  to  cases  wherein  U  is  infinite  dimensional  but  we  will  not  do 
that  here  . 

What  distinguishes  the  present  study  from  earlier  contributions  is  the 
attention  which  we  pay  not  only  to  the  relationship  between  the  operator.  A, 
and  die  input  element  b  ,  but  also  to  the  relationship  between  b  and 
the  semigroup  S(t)  generated  by  A  .  In  cases  where  A  has  discrete 


4 


spectrum  {X^,  j  k  e  K}  ,  K  being  a  countable  index  set,  this  amounts 
to  a  study  encompassing  the  input  element  b  ,  the  eigenvectors  {^j  k  e  K} 

of  A  ,  tv.e  corresponding  eigenvectors  of  the  dual  operator.  A'  as 

defined  in  Section  2  ,  and  the  exponential  functions  exp  ( X.^  t ),  k  e  K  . 

It  is  in  particular  reference  to  the  latter  that  what  is  probably  the  most 
important  idea  of  this  paper  is  developed  .  We  show  that  a  sufficient 
condition  for  b  to  be  an  "admissible  input  element"  ( definition  in  Section  2) 
can  be  given  in  terms  of  a  measure  on  Borel  subsets  of  the  complex  plane 
whose  support  is  e  K}  .  When  that  measure  turns  out  to  be  a 

Carleson  measure  the  input  element  b  is  admissible  .  This  result  brings 
out  yet  again  the  intimate  relationship  between  the  control  theory  of  infinite 
dimensional  linear  systems  and  parallel  developments  in  Hp  theory 
([5].  [8].  [12])  and  the  related  theory  of  completeness  and  independence 
of  sets  of  complex  exponentials  . 
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2.  Admissible  Input  Elements  . 

Let  X  be  a  separable  Hilbert  space  and  let  A  be  a  closed  operator 
on  X  with  domain,  g  (A)  ,  dense  in  X  ,  generating  a  strongly  continuous 

semigroup  of  bounded  operators  S(t)  on  X  for  t  >  0  .  For  b  e  X 

the  (generalized,  or  "mild")  solution  of 


x  =  Ax  +  bu ,  u  c  I*loc[0,  ®“)  , 
x(0)  =  xQ  e  X , 

is  given  by  the  M variation  of  parameters"  formula 

t 

x  (t)  =  S  (t)  xQ  +  f  S  (t  -  s)  bu  (s)  ds 


(2. 1) 
(2.2) 

(  3) 


and  may  be  seen  to  be  a  continuous  function  x  :  [  0,  ®)  X  .  Whether 
x  (t)  is  defined  for  each  t  2:  0  and  (2. 1)  holds  is  more  complicated: 
sufficient  conditions  are  that  b  e  £>(A)  or  that  u  is  differentiable  as 
a  function  of  t  ([3],  [11]). 

In  this  paper  we  wish  to  consider  (2. 1)  ,  (2. 2)  in  certain  cases  where 

b  does  not  lie  in  X  and  to  provide,  for  such  b,  a  formula  parallel  to 
(2.  3)  ,  Our  approach  is  similar  to  that  used  in  [14]  . 

I 

Identifying  X  with  its  dual  X  ,  we  denote  the  duality  relationship 
by  (x,  y)  ,  x  c  X ,  yc  X,  linear  in  both  x  and  y  .  Where  X  is 
the  complexification  of  a  real  Hilbert  space  X^  the  conjugate  element  "y 
is  well  defined  for  each  y  e  X  and,  with  (  ,  )  denoting  the  inner 

product  in  X  , 

(x, y )  =  <x,y>  ,  <x,y)  =  (x,y)  . 

The  bilinear  form  (  ,  )  is  symmetric,  i.  e. ,  (x,  y)  =  (y,  x),  x,y  e  X, 
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and,  for  all  x  e  X  , 

li x  II  x  = 


1  <x,y>  | 
sup  - 

r/o 


(2.4) 


The  symbol  A*  will  be  used  to  denote  the  dual  of  A  relative  to  the  bilinear 
form  (  ,  >  ,  that  is 


(Ax,  y)  =  (x,  A'y>  .  x  e  d(A)  ,  y  e  fl(A')  . 


The  operator  A'  is  closed  with  domain  jS>(A')  dense  in  X.  It  is  known 

that  if  A  generates  a  semigroup  S(t)  ,  then  S(t)'  is  also  a  semigroup, 
generated  by  A'  .  See  [  4  ]  for  details  . 


Let  Y  be  a  dense  subspace  of  X  which  is  a  Hilbert  space  in  its 
own  right  with  norm  |  |y  stronger  than  JJ  |JX  so  that  the  injection 
map 

J  =  Y  — ■  X 

j(y)  =  y  ,  y  e  y. 


is  one-to-one  and  continuous  with  dense  range  Y  C  X  .  We  further  suppose 
that  Y  is  invariant  under  the  action  of  S  (t)'  :  y  c  Y  =>  S  (t)'  y  €  Y  , 
and  that  this  map  is  continuous  with  respect  to  |J  S(t)*yUy  ,  (J  y  Uy  and  the 
usual  topology  of  [  0,  • )  . 

Let  Y*  be  the  dual  of  Y  with  respect  to  X  as  described,  e.  g.  in 
[1],  [W],  [15],  This  means  that  Y*  is  the  closure  of  X  with 


respect  to  the  norm 


sup 
y  €  Y 
y,l0 


|  <  x,  y  >  | 

II Y II Y 


II x  II  y*  = 


(2.5) 


It  is  known  that  Y  ,  so  defined,  is  a  realization  of  the  dual  space  of  Y 
and  it  is  easily  verified  that  the  bilinear  form  (  x,  y  )  may  be  defined,  by 
continuity,  for  x  e  Y4  ,  y  c  Y  as 

<x,  y>  =  lim  (x.  ,  y) 
k  —  ®  K 

where  {x^}  is  a  sequence  in  X  converging  to  x  in  ||  |y*  •  S 
defined,  ( x,  y)  generates,  as  x  ranges  over  Y*  ,  all  continuous 
linear  functionals  on  Y  ,  We  have 

X  C  X  C  Y1  . 

Definition  2. 1 .  In  the  system  (2. 1)  .  i.  e. , 

2 

x  =  Ax+bu,  u  e 

b  is  an  admissible  input  element  if  there  exist  Y.  Y1.  as  above,  with 


b  e  Y1,  such  that  for  every  T  >  0  the  continuous  map 


defined  by 


I*  :Y  -  C[0,T] 


(Ly)(t)  =  <b,S(t)'y>  ,  y  e  Y  ,  t  e  [0,T]  , 


(2.6) 


has  a  continuous  extension  to 


Lt  :  X  —•  LZ[0,T] 


Remark .  It  is  clear  that  this  amounts  to  the  statement  that  in  the  dual  observed 
system 


y  =  Ay 
z  =  <  b,  y  >  , 
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b  is  an  admissible  observation  element;  that  is,  for  yeY, 

z(-  )  =  <b,S(.  )'y  >  e  C  [  0,  T  ]  , 

this  relationship  extending  continuously  to  z  (•  )  e  L2[0,T]  for  y  e  X  . 

To  verify  that  Definition  2. 1  enables  consistent  definition,  at  least 
in  a  generalized  sense,  of  solutions  of  (2, 1),  (2.  2)  when  b  is  an 

admissible  input  element  and  to  establish  some  of  the  properties  of  the  resulting 
solution,  we  present 

Theorem  2.2.  If  b  is  an  admissible  input  element,  the  formula 

t 

(x  (t),  y)  =  <  xQ,  S(t)'y)  +  f  (b,  S(t-s)  y  >  u(s)ds  ,  yeY,  (2.7 

0 

defines,  for  each  t  a  0  ,  ji  unique  element  x  (t)  e  X  .  Given  T  >  0 
and  u  c  L2[0,T] 

x(t)  =  S(t)xQ  +B(t)u  ,  t  c  [ 0, T ]  ,  (2.8 

where  B(t)  is  the  strongly  continuous  family  of  bounded  operators 
B  (t)  :  L2  [  0,  T  ]  —  X  given  by 

<B(t)u,y)  =/<  b,  S  (t  -  s)'y  >  u(s)ds,  yeY.  (2.9 

Jo 

Proof.  From  (2.  8)  and  the  fact  that  Y  is  dense  in  X  it  is  clear  that 

x (t)  -  S(t)xQ  s  i  (t)  =  B(t)u 

where,  for  yeY, 

<£(t),y>  =  /  <b,  S(t-s)'y>  u(s)  ds  . 

0 

Let  x  e  X  and  let  {y^}  a  sequence  in  Y  converging  to  x  with 
respect  to  ||  *  •  Since  b  is  an  admissible  input  element  the 


corresponding  functions  defined  by 

=  <b,S(t-s)'yk> 

2  2 

converge  in  L  [  0,  T  ]  to  a  function  h  e  L  [  0,  T  ]  .  Defining 

<£(t),x)  =  /  h  (t-  s  )u(s)ds  , 

0 

we  see  that  for  t  €  [0,T] 

||<6(t).x>|*|h|  2  lu|  , 

L  l  0,  T]  L*[  0,  T  ] 

^UMllx||xi|u[[ 

A  Lz[o,T] 

since  (cf.  (2.  6),  (2. 10))  h  =  LjX  .  Hence  £  (t)  e  X1  =  X  .  1 

also  gives 

lUwlix*  UMIMi  2r  ,  . 

A  1  ir[  o,  t  ] 

showing  that  for  t  c  [0,T],  B  (t)  is  bounded  with 

1  b (t> ||  s  ||ltII  • 

To  establish  that  £  (t)  is  continuous  in  t  for  each  fixed 
u  c  L2[0,  T  ]  (and,  hence,  that  B(t)  is  strongly  continuous  in  t) 

A 

OststsT  and  form,  for  y  e  Y 

A 

<£  (t)  -  £  (t),y  >  =/  (b,  S  (t  -  s)'y>  u(s)  ds 

0 

-/  <  b,  S(t- s)' y  )  u(s)ds  =  (with  x  =  s  -  ("t  -  t ) ) 

0  t 

*  /  (b,  S(t-  x)y  )  u  (  x  +  (*t  -  t))  dx  -/  (b,  S(t-  s)'y>  u(s)ds 


+/  (  b,  S  (t  -  s  )'y  >  u  (s)  ds  = 

''o 
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A 

t  t-t 

=  J  (  b,  S(t-  s)'y>  (u(s  +  (t-t))  -  u(s))ds  +  /  <  b,  S(t  -  s)  y)u(s)  ds 
0  0 

*  1IM  liy||x(l|u(*  -  u||  +lul  2r  A  ). 

x  ir[  o,t]  lz[  o,  t  - 1] 

Since  Y  is  dense  in  X  and  since  for  fixed  u  c  L2[0,T]  we  have 

/s  lim  llull  2r  ^  .  =  °» 

t  - 1  ^  0 '  IT[  0,  t  - 1  ] 


linK  11  u(.  +(t-t))-uJl  ?  lim  |  u  ( •  +  (t -  t))  -  u|j  =0 

t  t  t  LZ[0.t]  t  i  t  I*  [ 0, t] 

We  conclude  that  for  fixed  u  e  L^[0, T]  ,  and  t,  t  as  described, 

^U^)-  IU£>  ~  ?  (t)llx  =  0 

t  t  t  A  t  *t 

and  thus  £  (t)  is  continuous  in  X  ,  This  completes  the  proof  of  the 
theorem  . 

Let  H  be  a  separable  Hilbert  space  and  let  {p^|  ke  K}  be  a 
sequence  in  H,  K  being  a  countable  ordered  index  set  .  The  p^  are 
strongly  independent  if  no  p^  lies  in  the  closed  span  of  {p^  J  1  /  k}  , 

If,  in  addition,  there  is  a  positive  number  c  such  that  whenever 


P  =  E  (2.11) 

*0 

the  being  complex  and  KQ  an  arbitrary  finite  subset  of  K  ,  we 

have 

E  |ak|2  a  o2|lp|l|  (2.1 

Ko 

2 

we  say  that  the  pfc  are  uniformly  f  -independent,  since  (2,12)  implies 

E  l«kl2  *  °2IpIh 

K 


(2. 12) 


(2. 13) 


11 


whenever  { cl2  and  p  =  ^  is  convergent  in  H  . 

K 

If  there  is  a  positive  number  C  such  that 


pIh  3  °2  E  ki: 


p  as  in  (2.11),  we  say  that  the  sequence  {p^}  is  uniformly  f  -convergen 


Since  this  property  implies  that  if  {a^}  6  *2  the  series  E<y]cpic  ls 

K 

convergent  in  H  and 


I! Pi!  a  c2£  |ck|2. 


(2. 14) 


Recall  that  a  sequence  {  p^}  in  H  forms  a  Schauder  basis  for 
H  if  for  every  p  €  H  there  are  unique  coefficients  such  that  the 

series  Za]cp)c  converges  to  p  in  H  ([21]).  A  Schauder  basis 


which  is,  at  the  same  time,  both  uniformly  I  -independent  and  uniformly 


I  -convergent  is  a  Riesz  basis  .  For  evident  reasons  we  shall  also  use. 


2  2 

synonymously,  the  term  uniform  I  -basis.  If  {p^}  is  a  uniform  l  -basis 


for  H  then  every  p  in  H  has  a  unique  convergent  representation 


*  £  “i 


with  (cf.  (2. 13),  (2.14)) 


<='2E|°vIz  =  Up||: 


c2£  i«ki2. 


For  the  remainder  of  this  section  we  suppose  that 


(1)  the  operator  A  with  dense  domain  fi(A)  c  X  generates  the  strongly 


continuous  semigroup  of  bounded  operators  S(t),  t  »  0  ; 


consists  of  discrete,  simple 


(ii)  a  (A)  ,  the  spectrum  of  A  , 

eigenvalues  X^,  k  e  K  ,  and  the  corresponding  normalized  eigen¬ 
vectors  k  e  K  ,  form  a  strongly  independent,  uniformly 

convergent  Schauder  basis  for  X  . 

Since  the  fi  ,  k  €  K  ,  are  strongly  independent  and  have  closed 
span  equal  to  X  there  exist  unique  biorthogonal  elements  i(j^,  k  e  K, 


such  that 


f  1,  *  =  1 

^k*  =  \  ,  k,  I  e  K  . 

k  1  \o,  M  1  , 

As  is  well  known,  the  ^  are  eigenvectors  of  the  dual  operator  A' 
corresponding  to  the  eigenvalues  X^,  k  e  K  .  We  further  assume 
(ill)  the  eigenvectors  of  A'  have  the  property 

\  €  Y  C  X  . 

(this  is  true,  for  example,  if  YD  S((A')r)  for  some  positive  integer  r) 

If  x  e  X  ,  the  fact  that  the  ^  ^  form  a  Schauder  basis  in  X 

implies  the  existence  of  unique  £ ^ ,  k  e  K  ,  such  that 

K 

the  series  converging  in  X  ,  From  this  it  is  evident  that 

=  (  vhc»  x>  »  k  €  K  . 

We  are  not  assured,  in  general,  that  the  are  square  summable  but  the 

uniform  i  -convergence  property  of  the  ^  shows  the  square  summability 
of  the  sequence  to  be  a  sufficient  condition  for  convergence  of  (2. 15)  . 


Since  we  assume  the  4^  lie  in  Y  ,  given  any  element  b  c  Y1 


( and  this  includes  b  e  X )  we  may  define 


bk  =  (  ^  ,  b  >  (2. 16) 

and  obtain  a  set  of  coefficients  ,  ke  K,  associated  with  b  .  In 
general  it  is  not  possible  to  recover  b  from  the  coefficients  .  (An  examp 
is  X=L2[0,  2ir]  ,  Y  =  H^[  0,  2i]  ,  i^fx)  =  (2ir  fV**  ,  k  =  0, 1,  2,  .  .  .  . 
The  ipjJ  *  ^k)  here  form  an  orthonormal  basis  for  X  and  belong  to  Y 


but  there  is  a  non-zero  element,  namely  6/rt.  -  6.,  .  , 

(0)  (zir) 


in  Y  for  which 


all  of  the  b^  are  zero.  This  arises,  of  course,  because  the  closed  span 

of  the  in  Y  is  not  equal  to  Y  . )  As  a  consequence  it  is  not 

generally  meaningful  to  write  b  *  E  • 

X 

Nevertheless  it  may  be  meaningful  to  consider  the  initial  value  problem 
(2.1),  (2.2),  i.  e. , 

x  =  Ax  +  bu 

x(0)  =  x0£X,  ueL2toc[0,-), 

for  certain  b  c  Y'  ,  namely,  those  that  we  have  already  characterized  as 
admissible  input  elements  .  We  wish  now  to  show  that  the  class  of  such 
admissible  input  elements  can  be  characterized  in  terms  of  the  coefficients  b^ 
and  the  eigenvalues  X  k  .  If  x  (t)  is  the  solution  of  (2. 1),  (2. 2) 
established  by  Theorem  2.  2  for  an  admissible  input  element  b  then, 
in  particular,  for  t  a  0  , 


(xft),^)  =  <xQ,  Sft)'^) 


+  /  (b,  S(t-s)‘  4^}  u(s)ds 

V°  Vt"s) 

=  e  x0Jc+J  e  u  (s)ds 
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where 


x0  =  E  xo,  k  ^  k  * 
K 


We  do  not  know  that  the  numbers 
but  the  series 


V 

e  x, 


0,  k 


are  square  summable 


K 

must  converge  to  S(t)xQ  by  virtue  of  the  (assumed)  Schauder  basis  property 
of  the  .  It  follows  that  a  sufficient  condition  for  x  (t)  to  belong  to  X 
is  that  the  numbers 

rt  X.  (t-s) 

Sk(t)  =  bkjT  e  u(s)ds  (2.17) 

should  be  square  summable  for  each  t>0.  Equivalently,  making  a  trivial 
change  of  independent  variable , 

, t  V.s 

=  e  f(s)ds»  f(s)  =  u(t-s)  . 

The  necessity  of  considering  an  infinite  number  of  values  of  t  can  be  obviated 

by  taking  f  to  be  an  element  of  L2[0,T],  T>0  fixed,  and  defining 

f  (s)  s  0  in  [t,  T]  for  t  <  T  .  The  map 

T  X  s 

Sk  =  \  6  k  f(s>ds»  f  c  L2[0,T], 

so  defined  may  be  designated  as 

li  *  l2[o,t]  -  X, 

LT(f)  =  x  =  E  Ck  ^k  ' 


(2. 18) 

(2. 19) 
(2.  20) 


and  it  is  easy  to  see  that  L't  Is  the  dual  of  :  X  — ■  L2  [  0,  T  ]  as 

defined  by  (2.  6)  .  Thus  the  boundednes  s  of  L,^  , 


as  required  in 


Definition  2. 1,  may  be  obtained  as  an  immediate  corollary  if  it  is  shown  that 
L't  ,  defined  by  (2.18)  -  (2.19)  ,  is  bounded.  For  our  present  purpose 
this  is  the  route  of  choice  . 

Extending  f  further  via  f  (t)  =  0  ,  t  >  T  ,  the  Laplace  transform 

of  f  is  the  entire  function 

®  T 

t(z)  =  /  e~zt  f(t)dt  =  /  e"ztf(t)dt. 

0  0 

In  terms  of  we  clearly  have 

Ck-M  K- 


and  the  following  proposition  is  evident  . 


E?l!ls!aSSi5a 


L^,  )  is  bounded 

lust  in  case,  for  every  fe  T]  the  Laplace  transform  of  f ,  ip  has 

the  property 

Elbk^("\)|2  <  »•  (2.21) 

K 

We  are  fortunate  that  the  inequality  can  often  be  established  with  the 
use  of  the  concept  of  a  Carleson  measure  and  the  corresponding  Carleson  measure 


theorem  as  it  applies  to  the  space 


=  H^{z  |  Re(z)  >  a)  ,  a  real  . 


(2.  22) 


2 

The  space  H  {z  |  Re(z)  >  a}  consists  of  those  complex  functions  ^( z), 
analytic  in  Re(z)  >  a  ,  bounded  in  each  half  plane  Re(z)  >:  a  +  6,  6  >  0 
and  satisfying 


/  +m)l2dn  *  M,  ,  $  >  or  , 


(2.  23) 


where  is  a  positive  number  depending  only  on  ( and  not,  in 

particular,  on  £  )  .  It  is  known  (  see,  e.  g.  [  10]  )  that  each  such 
function  has  a  limiting  "boundary"  function 

J^)  =  Mm  ^  +  iH)  (2.24) 

£  *  a 

defined  almost  everywhere  in  -»  <  T)  <  <»  and  <b  (t))  is  measurable 

•  a 

with 

f“  l^ao)|2  <in  s  . 

2  2 
Each  €  Hfl  is  the  Laplace  transform  of  a  unique  function  f  e  L ^  [  0,  «  ) 

such  that 

/  |e~atf(t)j2dt  <  .  . 

0 

Let  p.  be  a  (non-negative  valued)  measure  defined  on  the  Borel 
subsets  of  {z  |  z  >  a}  .  Then  p.  is  a  Carleson  measure  if  for  every  real 
t  and  every  h  >  0 

p.({z|T-h5lm(z)s  x  +h,  a  <Re(z)  s  a  +  h»  sAh  (2.  25) 

for  some  positive  A  depending  only  on  p.  ( not  on  h )  . 

For  a  Carleson  measure  we  have 

Theorem  2. 4 .  If  p.  is  a  Carleson  measure  on  {z  j  Re(z)  >  a}  with  A 
55.  in.  (2.  25)  ,  if.  H2  ,  and  is  given  by  (2.24),  then 

f  U(2)|2dH(z>  3  i220A  /”  |jrf  (n,|2  dn  .  (2.26) 

{z|Re(z)>  at}  it  -« 
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A  proof  of  this  theorem  is  offered,  for  the  sake  of  completeness,  in 
Section  4  of  this  paper  .  The  relevance  of  this  theorem  for  our  present  studies 
is  exhibited  in  the  selection  of  a  particular  measure  p.  .  For  be  Y'  and  a 
given  discrete  spectrum  {X.^}  for  A,  let 

1  11  * 
be  defined  by 

=  Ibk|2  ,  keK,  (2.27) 

|  p({z|Re(z)  >  a}  -  {\  |k  e  K})  =  0  .  (2.28) 

In  this  case  the  left  hand  side  of  (2.  26)  becomes 

1 

Z  i*>k**<-xkH2  <cf*  <2-2l»  • 

K 

The  Plancherel  Theorem,  on  the  other  hand,  gives 

/  If*  (H)|2  dTl  =  2u  /  |e_atf(t)|2dt 
-<■  0 

^  2e2MTir  fT|  f  (t) |2 dt 
0 

when  the  support  of  f  is  restricted  to  [  0,  T]  .  Thus 

Z  200°e2|“|T  $  /T  U (t)|  2  dt 

K  0 

and,  in  view  of  our  earlier  discussion,  we  have 

Csa2yj|iy=|±_| .  A  sufficient  condition  in  order  that  b  e  y'  should  be  an 

admissible  input  element  for  the  system  (2.1)  ,  wherein  a(A)  =  {X^  |k  e  K} 

and  the  corresponding  eigenvectors  k  e  K,  form  a  strongly  independent  , 

2 

uniformly  C  -convergent  Schauder  basis  for  X  ,  is.  that  the  measure 
U.  «■»  i  defined  by  (2.27),  (2.28)  should  be  a  Carleson  measure  in 


v  •,  '  *«  *«, 


{  z  |  Re  (z)  >  a)  for  some  real  a  . 

We  remark  that  the  assumption  (i)  above  together  with  the  Hille-Yoshida 
Theorem  ([4],  [11])  implies  that  the  complex  numbers  -X^,  ke  K,  are, 
indeed,  confined  to  some  right  half  plane  Re(z)  >  a  .  The  fact  that  the 
support  of  f  is  restricted  to  [0,T]  implies  that  the  corresponding  Laplace 
transform  ^  is  entire  and  satisfies  an  inequality  (2.  23)  for  every  real  a 


Identification  of  Admissible  and  Inadmissible  Input  Elements:  Examples 


Our  first  task  in  this  section  will  be  to  develop  a  method  whereby  input 
elements  b  not  in  the  state  space  X  may  be  identified  as  particular  elements 
of  a  larger  space  Y1  .  The  assumptions  made  will  be  somewhat  more 
restrictive  than  thos  introduced  in  Section  2  .  They  are  by  no  means  necessary 


conditions  . 


Let  us  suppose  that  the  operator  A 


generating  a  strongly  continuous 


semigroup  S(t)  on  the  Hilbert  space  X ,  has  (dense)  domain  fl(A)  and 


that  A  possesses  discrete  eigenvalues 


k^ ,  k  e  K  ,  with 


lim  Ik.  I  =  ». 
pOc)-«  K 

Here  p  (k)  denotes  the  number  of  elements  I  e  K  such  that  I  <  k  with 
respect  to  the  assumed  order  relation  on  K  .  The  corresponding  normalized 
eigenvectors  ^  are  assumed  to  form  a  uniform  basis  for  X  .  We  denote 
the  dual  operator  by  A*  .  It  has  the  same  eigenvalues  k^  and  the  correspond! 
eigenvectors  4^  »  k  e  K  ,  will  be  assumed  normalized  so  that 


<V*I>  ■ 


1  ,  k  =  f 

0,  k  ^  f  . 


The  4^  also  form  a  uniform  basis  for  X  ,  as  is  well  known  .  Then  it  is 
easy  to  see  that 

fl(A)  =  {y  =  Z  xk  Z  I  i2  <  “  * 
s  K  K 

and  that 

fl(A')  =  {y  =  ^y  4»klZ  I xkykl2  <  **  • 

K  K 

For  the  work  of  this  section  we  take  Y  =  jB(A')  with  the  graph  norm 


20 


where 


UyIIy  s  Ea  +  lxicl2)  |ykl2* 

K 

OB 

Y  =  E  Yjc^j, 


k  =  1 


in  X  .  Then  Y  C  X  and  the  injection  mapping  is  continuous  .  It  will 
often  be  possible  to  identify  a  Hilbert  space  Z  C  X  with  continuous  injection 
map  such  that  j]  J|^  is  a  familiar  (e.g.  Sobolev)  norm  and  Y  is  a 

closed  subspace  of  Z  on  which  the  norms  ||  §  z  and  ||  JJ  y  are 
equivalent  . 

We  will  be  concerned  with  two  different  extensions  of  the  operator  A  . 

A 

We  suppose  first  of  all  that  there  is  an  element  x  e  X  not  in  fl(A)  and 
that  L  is  an  operator  on  X  such  that 

6(1)  =  {£  +ux  |  £  c  fl(A),  u  scalar}  , 

Lx  =  Ax,  x  e  j9(A)  . 

We  will  refer  to  L  as  an  "operational  extension"  of  A  .  Its  significance 
arises  from  the  fact  that  many  of  the  inhomogeneous  boundary  value  problems 


arising  in  applications  can  be  expressed  in  the  form 

dx  _ 

dt  ' 


(3. 1) 


with  the  restriction 


x  =  5+uxe  «(L)  . 


(3.  2) 


The  second  extension  of  A  ,  which  is  a  map 

£  :  X  —  Y’  , 

is  a  standard  one,  often  used,  e.  g.  in  [14]  .  If  y,  T  e  fl(A)  ,  $(A')  , 

respectively,  we  have 


<  Ay,  n  )  =  <  y,  a'  n  >  . 

Since  A'  :  Y  =  fl(A*)  —  X  is  continuous,  the  form  (y,A‘  H>  extends 
to  (x,  A'n)  .  xcX,  by  continuity  and  density  of  fi(A)  in  X  and, 
so  extended,  ( x.  A1  *1  >  defines,  for  each  fixed  x  e  X  ,  a  continuous 
linear  functional  on  Y  ,  i.  e. ,  an  element  of  Y1  .  We  define 

A  :  X  — ■  Y1  =  (  A(A'))' 


(Ax,  T]  >  =<x,A,n),  xeX,  Ti€Y=a(A'). 

Our  first  goal,  with  reference  to  the  system  (3.1),  (3.2),  is  to  replace 


it  by  an  infinite  set  of  scalar  ordinary  differential  equations 

ET  *  Xk*k  +  b]cU'  k  e  K’ 

where 

x(t)  =  ^  XjJt)  , 
ke  K 

convergent  in  X  .  In  order  to  do  this  we  recognize  first  of  all  that 

z  -  2  WK 

K 

A 

represents  not  Lx  ,  but  rather  Ax  ,  since 

=  (X,  a\)  =  <x,  =  ^.x,^  . 

We  rewrite  (3. 1)  in  the  form 

~  =  Ax  +Lx  -  , 

at 

an  equation  in  y'  .  Then,  since  x  is  to  have  the  form  (3. 2)  with 
£  e  £(A)  ,  and  since 


(3.3) 


(3. 4) 


L£  =  A£  =  A|  ,  §  €  fi(A) 


(3.4)  becomes 

Hv  /V  A  AA 

^  =  Ax  +  (Lx  -  Ax)  u  . 

We  define  b  c  Y*,  a  continuous  linear  functional  on  Y  =  fl(A'),  by 

<b,n>  =  (i£  -  ax,  n  >  =  <l£,  n  )  -  (x,  a'  n>  (3.5) 

for  t]  c  fl(A')  =  Y  .  We  then  have 

b=  EMk 

K 

where  the  "control  input  coefficients",  bk  ,  are  given  by 

bk  =  <b,4jk>  =  (l£,  ijjk>  -  (x.A'^) 

=  (lx,  4^)  -  x-jc(  +k>  •  *3*6^ 

In  most  examples  we  shall  have  Ijc  =  0  .  Then,  if 

K 

convergent  in  X,  we  obtain,  in  place  of  (3.6)  , 

bk  =  "Xk  *k  *  k  e  K  •  (3*71 

Also,  in  this  case,  the  equation  (3.  5)  becomes 

<b,  T)  )  =  .(  x.  A*  H>  .  (3.  8] 

The  equation  (3. 5)  (or  (3. 8))  will  generally  be  used  to  identify  the 
functional  form  of  b  while  (3.6)  (or  (3.7))  will  be  used  to  identify  its 
expansion  coefficients  in  terms  of  the  eigenvectors  ^k  of  the  operator  A  . 
While  not  all  admissible  input  elements  can  be  treated  this  way  the  class  is 
large  enough,  we  believe  to  warrant  the  detailed  description  we  have  given  here. 


Example  1 :  Heat  Equation .  Let  x(s,t)  satisfy 


9x  a  x 
at  ~  “772 


o  <  s  <  l ,  t  >  o  , 


with  boundary  conditions 


x(0,t)  a  0  ,  “x(l,t)tP|.(l,t)  =  u(t)) 
where  a,P  are  real  numbers,  not  both  equal  to  zero  .  In  this  case 


we  take 


Ax  = 


LX  a 


x-  I*  i o»  1  ] 

x  c  *(A)  =  {x  e  H2[0, 1]  j  x(0)  =  0  , 
ax  (1)  +  Px'  (1)  =  0}  , 

x  €  fi(L)  =  {x  c  H2[0,1]  |  x(0)  a  0 >  , 


A 

X  (s)  = 


With 


’■oTBS-  a+p  *  0 

,  a  +  P  =  0  . 


<x,y  >  =  f  x(s)  y  (s)  ds 
0 


we  see  that  if  x,  y  €  S(A) 


(3.9) 


(3. 10) 


(3. 11) 


<Ax,  y>  -  <x,Ay)  =  f  (x"(s)y  (s )  -  x  ( s  )  y  "(s))  ds 

0 

d  t 

ds  (x  (s)  y(s)  -  x(s)y'(s))  ds  =  (since  x(0)  =  y(0)  =  0) 
x'(l)y(l)  -  x(l)y'(l) 

/(x'(l)  +  |x(l))y(l)  -  x(l)(|y(l)  +y'(l)),  0^0 

\x’(l)(y(l)  +  1/(1))  -  ( x (1)  +!■  x'(l))y'(l) ,  a  {  o 


and  we  conclude  A  =  A  .  In  the  first  case  of  (3. 11)  ,  a  +  P  ^  o  , 

A  t 

Lx  a  0  and  we  have,  for  n  e  £(A  )  =  £  (A) 

<b,  n>  =  -  <x,  A'n)  =  J  s  n" (s)  ds 


!  ,i  i  -n’(i)  +  T)(i) 

■  aT@  (-s1(s>l0  +,><s)l0) - o-Tfj— 


Thus  we  have 


~  na).  Mo 

p 

-S-n'(l),  a  /  0  . 


I6(l)'  0 

1  6'(1)  ,  a  ^  0  . 


The  two  agree  if  neither  c i  nor  P  are  zero  because  the  linear  functional 

^  5(1)+  S’  is  zero  in  (  S(A))'  =  Y1  in  this  case  . 

2 

The  eigenvalues  of  A  are  where,  for  k  =  1,  2,  3,  . . 

a  sin  (u>.)  +  Pco.  cos(«.  )  =  0  .  j 


=  sin  0, 


and  (3.13)  becomes 


so  that 


'  2  ,o2  2 

a  +P  01 


=  cos  6 . 


cos  (u^-  ek)  =  o 


g<icafc 


SS 

:« 


£.3 


K 


•v 


;*v 

£5 


M 


rs 

$ 

Vt 


$ 

b 


Si 


,2k-  1  ,  .  ,  -1 

oj  .  =  ( — = — )  it  +  sin 


a 


Jn 2  .  ft2  2 
Wa  +  p  a> . 


It  is  easy  to  see  that  — ®  (  i-)  as  k  — 

wk  k 


so 


01, 


=  (-%^-  )  ir  +  {>  (^  )  , 


,  2k  -  1  ,  .it  , 
“k  (— y~  +  2  s  kir 


P  /  o, 

P  *  0. 


(3. 1- 
(3.1! 


Defining 


vk=  £  sln<“ks)' 


ds 


it  is  easily  seen  that  in  all  cases  the  vk  are  nonzero  and 

Hm  v  =  -A-  . 

k  — «  K  >JT 

Then  the  eigenfunctions 

f*k(s)  =  ~v~  sin(o>ks) 

2 

form  an  orthonormal  basis  for  L  [0, 1]  .  It  follows  then  that  the  coefficients 
of  the  input  distribution  elements  (3.12)  are  given  by 


bk  ■ 


Pv, 


sin 


<“k> 


0  it  0. 


«k  COS  («k) 


(3. 


a  v. 


a  ^  0  . 


We  consider  here  the  case  P  -ft  0  ,  saving  the  analysis  for  P  =  0 
until  later  in  this  section  .  If  P  /  0,  formula  (3. 16)  shows  the  bk 

to  be  uniformly  bounded  .  The  complex  numbers  -  ^  =  w k  have  the 

property  (from  (3.14)) 

~\  =  ( — 2 — *  *  +  * 


(3. 


Thus  the  number  of  such  in  any  set  -  j  Im  (z)  -  t|  s  h  , 

Re  (z)  2S  a  +h  is  ®(h^^)  and  it  follows  that  the  measure  n  with 

CD 

H(-\>  =  \\r  .  t*({Re(z)  -  a}  -  U  )  =  0  is  a  Carleson  measure. 

Hence  if  fr  0  the  boundary  input  (3.10)  is  admissible  . 

In  this  case  the  result  is  easily  obtained  without  the  Carleson  measure 
theorem;  for,  if  the  coefficients  c^  are  square  summable  and  T  >  0  , 


*  v 

IE  bkclce^llI2rn  T1 

k  =  l  I*  [0.T] 


supdbfcl}/  E  I  <^1' 

*  k  =  1 


J 


v  rT  2  V 
E  /  e  *  dt 

k=  1  0 


sup 

k 


1‘tWe  /*e2Xfct  dt  =  sup{|b|>  J  1  l^l2 

k- 1  k=  1  0  k  k=  1  k=l 


-1 


Nc 


(3.18) 


since  sup{|bjcJ}< 


and  we  conclude  that  the  function  sequence 


.-v. 


{b^e  *  }  is  i  -convergent  in  L[0,T]  .  Our  next  example  is  chosen 
in  such  a  way  that  a  simple  argument  of  this  type  does  not  apply  and  the  Carleson 
theorem  is  actually  needed  . 


Example  2.  Another  Heat  Conduction  System  .  As  a  further  example 


we  ask  the  reader  to  consider  the  system  shown  in  Figure  3. 1 


Figure  3. 1 


The  shaded  horizontal  bar,  B  ,  represents  a  layer  of  material,  whose  depth 
will  be  assumed  negligible,  and  whose  heat  conductivity,  k  ,  is  small  in 
comparison  to  its  specific  heat  R  while  the  region  n  consisting  of  the 
half  strip 

fi  s  0  ^  x  :s  1,  z  ss  0  , 

is  assumed  filled  with  a  material  whose  specific  heat,  r  ,  is  small  by 
comparison  with  its  conductivity,  K  .  The  heat  flow  equations  are  thus 
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R  3T  _  .  K  a  t 

*  ax2  az  ’ 


(3. 19) 


ix  „  .a2?  .  a2x 
at  “  K  { 


2  +  2 

ax  az4 


)  , 


(3.  20) 


together  with  boundary  conditions 


4|<o.t)  =  o,  H<i.t)  =  o 

■§£  <0.z,t)  =  0  ,  =  g(z)u(t). 


11m 

Z  -*  -« 


~(x,  z,t)  =  Urn  r(x,z,t)  =  o 

“  Z  — •  —OB 

x(x,0,t)»  T(x,t),  Osxsl  . 


(3.  21) 
(3.  22) 
(3.  23) 
(3.  24) 


The  Inhomogeneous  boundary  condition  along  x  =  1,  z  ^  0  ,  represents 
the  input  heat  flux  .  In  (3.19),  (3.21),  T(x,t)  is  the  temperature  in 

the  bar,  t(x, z,t)  the  temperature  in  O  . 

If  we  assume  k,  r  vel'y  s;.idU  by  comparison  with  R,  K,  we  may, 
as  an  Idealization,  replace  f  3-.  19J  and  (J.  20)  by 


p  8t 

at  ” 


82t 


+  -^4-=  0, 


2 

ax**  8z‘ 

retaining  the  boundary  conditions  (3.21)  -  (3. 24)  . 
state  space 

X  *  {T  *  T  (x)  j  T  e  L2[0,1]  >  . 

We  define  an  operator  A  on  X  with  domain 

8(A)  a*  H1  [0»  1] 


(3. 25) 


(3. 26) 


We  take  as  our  basic 
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as  follows  .  Given  Tc  A  (A),  we  let  t  =  t(x,  z)  satisfy  (3.26) 
in  ft  together  with 


t(x,0)  s  T  (x),  05XS1, 

(3. 27) 

and 

|i.(0,z)=  0,  -fj  d,z)  »  0 

(3. 28) 

lim  «|f  (  •  ,  z)  a  0  in  L2[0,1]  . 

2  — » 

(3.  29) 

lim  T(*  ,  z)  a  0  in  H^O.l]  . 

Z  —  -oo 

(3.  30) 

From  [14]  , 

for  T  €  H*[0, 1]  we  have  re  H3^2(ft)  .  The  trace 

theorem  ([  1  ] 

,  [14])  then  gives 

^f  (.,0)€  L2[0,1], 

and  we  define 

AT  =  -  |  |f  (•  ,  0)  .  (3.31) 

So  doing.  (3.25)  becomes 

T  a  AT  (3.32) 


and  (3. 31)  is  subsumed  In  the  definition  of  A  . 


Lemma  3.1. 
operator 


The  operator  -A  is  the  positive  square  root  of  the  Sturm- Liouville 


with 


•s- 


i(T)  =  {T  '€  H  [  0,  1]  |  Hi(0)  =  (1)  =  0}  . 

Proof.  We  compute  (-A)2T  for  T  e  fi(T)  .  Far  such  T  the  solution 
of  (3. 27)  -  (3.  30)  €H5/2(0).  If  we  let 


then 


~<x’z>  =  I  Ti  <x’z> 


T  (  • ,  0  )  a  -  AT 


( -A)2  T  =  -At  (•  ,  0)  =  ^  (•  ,  0) 

IT  8z* 

-  K2  dfr  .  _  K2  d2T 

R  8x  K  dx 

since  t€  H5/2(SJ)  together  with  (3.26)  implies  that 

2  2 

^-4(-  ,  0)  +  -^(  -,0)  =  0  in  L2[0,1] 

ax  az* 

and  T  *  t (  •  , 0  )  . 

The  positivity  of  -A  follows  from  the  divergence  theorem  .  If 
T  €  fi(A)  and  t  =  t(x,  z  )  is  constructed  as  above,  we  have 

If  1<&<X-Z»Z  +  (4j(x,z))2]  <*dz 


=  //  JVt(x,  z  )1(2  dx  dz  (V=  gradient) 

n 

a  //[div( T(x,z)Vx(x,z))-t(x,z)A2T(x,y)]  dx  dz 


(AT  -  Laplacian)  =  (from  (3.26)) 

a  ff  div(x(x,  z )  Vt(x,  z ))  dx  dz  a  (using  (3.  27)  -  (3.  30)) 

*  f  T(x,0)|?  (x,0)dx=  (T,  -AT)  , 


I 


I 


This  completes  the  proof  . 

Accordingly,  A  is  self  adjoint  with  eigenfunctions 

fi  q(x)  si,  ^  (x)  =  sfz  cos(kirx),  k  =  1,2,3,... 
and  eigenvalues 


x0  =  °*  Sc  =  “  TT  *  k=l,  2,  3, .... 

Let  w  (x,  z )  be  the  solution  of  the  following  inhomogeneous 
boundary  value  problem: 


82w  ,  a2w  »  , 

- _  +  - -  =o  in  Cl 


8xfc 


8z 


I^(0.z)«0,  ( 1,  z )  =  g (z) 

lim  -|^-(x,z)»  11m  w(x,  z)  =  0  . 

Z  —  -OD  z 


w  (x,  0 )  s  0  , 


1 . 


We  will  assume  that  g  (z)  is  such  that  the  resulting  w  (x,  z)  t  H' 


In  this  case  the  inhomogeneous  equation  can  be  interpreted  as 
T  =  AT  +bu 


where  b*  b  (x)  is  given  by 

b(x)  s  “I  I F*(x’°>  •• 

To  compute  the  coefficients  of  the  expansion 

09 

b  (x)  =  Yj  \  **k  *x)  * 
k  =0 


we  note  that  since  A  is  self  adjoint,  4^  (x) 

bks  f  **k*x*  b  M  **  • 


and 


i*yy. 


Let  $  c(  x,  z )  be  the  solution  of 


d2*y.  a2^ 

- +  - =£-  =0  in  « 

ax'2  az^ 


with 


<&k(x,  o)  a  <&k(x) 


and  homogeneous  boundary  conditions  of  the  type  (3.  27)  -  (3.  30)  otherwise. 

Then,  with  a2  =  +  -K-  , 

8x*  8z* 

r  2  2 
o  =  f  [<tk(x,  z  )  A  w  (x,  z)  -  w  (x,  z)  a  ^  (X,  z )]  dx  dz 

=  f  div[$k(x,z)  grad  w(x,z)  -  w(x,z)  grad  ^(x,  z )]  dx  dz 


R  r1 

"  "K  f  4MbMdx  +  /  ®1C(1»Z)  9(z>  dz 


giving  (cf.  (3.35)) 


K  r° 

hks  R  f  *k  ( 1.  z  )  g  (z)  dz  . 


Now  it  is  easily  checked  that  for  k  =  1,  2,  3, . . . 


so  that 


and  thus 


^(x,  z )  s*  (>JT  cos  kirx)  (exp  (kirz )) 


<&v  (1.  z )  =  (-1  r  */T  exp  (kir  z ) 


fcL  =  K  J  exp{kirz) 
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The  Carleson  measure  theorem  can  be  used  in  a  slightly  different  way 


than  that  set  forth  in  Corollary  2.  5  to  show  that  if  g  e  L  ( -® , 0 ]  then  the 


2 

d,.  are  square  summable  and  b  is,  consequently  an  element  of  L  [0,1]. 


Writing  £  =  -z  ,  g(-£)  =  g(£),  we  see  that 


b  =  f  exp(-kw£)  g(£)d£  . 


Since  the  measure  p.  assigning  the  value  1  to  each  of  the  points  ktr  , 

f-l)k^  K 


k  =  0,1,2,.,,,  is  clearly  a  Carleson  measure,  and  since 

.2 


changes  only  in  sign,  {b^}  e  l‘ 

If  g  (z)  is  just  bounded  and  measurable  on  -®  <  z  s  0  we  can 
llmost  trivially  obtain 

bfc  =  S(^) 

and  the  b^  will  be  square  summable  . 

It  is  obviously  possible  to  replace  g  (z)  by  distributions  of  various 


types  .  Taking  g  (z)  =  6 


(0) 


corresponds  to  a  point  heat  source  at  the 


comer 


x  =  1 ,  z  =  0  and  leads  to 


.  =  j-C.g.K 

°k  R 


(3.  36] 


2  1 

In  our  present  example  X  =  L  [  0, 1]  ,  Y  =  £  (A)  =  H  [0, 1]  and 


r-1 


Y  =  H  [0,1]  .  The  coefficients  (3.36)  may  be  recognized  as  those 


corresponding  to  6 


(1) 


( referring  now  to  distributions  along  the  x-axis)  . 


Any  measure  p.  assigning  to  the  points  ^  ktr  values 


jbjJ  which  are  bounded  evidently  yields  a  Carleson  measure  and  we  conclude 
that  all  of  the  above  cases  correspond  to  admissible  input  elements  .  In  this 
case  the  argument  represented  by  the  inequalities  (3. 18)  will  not  work 


because  the  series 


I  c- 


-1 


TV 


)  is  not  summable  in  this  example 


Example  3 


Hyperbolic  and  Neutral  Systems 


A  wide  variety  of  systems  involving  linear  hyperbolic  partial  differential 
equations  in  two  independent  variables  x,  t,  or  neutral  functional  equations 
lead  to  systems  of  the  form  described  at  the  beginning  of  this  section,  the 
eigenvectors,  ^  ,  of  A  forming  a  uniform  i2-basis  for  the  state 
space  X  and  the  eigenvalues  confined  to  a  vertical  strip 

a  <  Re  (X )  <  P  in  the  complex  plane.  It  also  usually  turns  out  in  these 
cases  that  the  number  of  X^  in  any  rectangle 

a  <  Re  (X )  <  P,  y  <  Im(X)  <  6 

Is  less  than  or  equal  to  M(6  -  y)  ,  where  M  is  a  fixed  positive  number. 
It  is  evident  that  the  measure  (2.  27),  (2.  28)  is  a  Carleson  measure  in  tfo-ae 
cases  whenever  the  control  input  coefficients  b^  constitute  a  bounded  set  . 

Example  4  .  Linear  Surface  Waves  .  If  the  operator  A  is  defined  as 
in  (3.  31)  but.  Instead  of  the  first  order  system  (3.  32)  we  consider  the  second 
order  counterpart 

£  +  A  £  =:  o  (3.  3 

we  obtain  the  linearized  equations  for  small  amplitude  waves  on  the  surface 
of  an  incompressible  fluid  .  The  theory  is  more  fully  developed  in  [16]  ,  [17]  , 
[19]  .  With  *1  =  £  ,  (3.  37)  is  equivalent  to  the  first  order  system 

(n)=  (-a  o)(n)s  G(n)  •  <3-3 

To  obtain  a  topology  corresponding  to  the  energy  of  the  system  one  defines 


where 


i\n  J» 


=  l?i2r  +  (n.A_in)  2r 

e  L 2  [0,1 ]  L2q[0,1] 


V- 


2  2 

LqCo.IJ  =  {<;  e  LZ[0,1]|/  £  (x)  dx  =  0 }  . 

0 

2 

The  restriction  to  LQ  [  0, 1]  corresponds  to  conservation  of  fluid  volume 
On  the  domain 

«0(A )  =  (C  e  H*[  0,  1]  |  ^  £  (x)  dx  *  0  } 
the  operator  A  Is  Invertible  .  Its  eigenvalues  are  (cf.  (3,  34)) 

Nc  =  “  R  ^ 11  ’  k  =  1,2,  3,,,, 

with  the  same  eigenfunctions  ^^(x)  ,  k  =  1,  2,  3, . . .  ,  as  shown  In 
(3. 33)  .  Correspondingly,  the  operator  G  has  eigenvalues 


u>k.  -1<V  «>k  = 


-  Kir  1/2  vl/2  =  kl/2 


k'  “  =  ykA/0,  k  =  1,2,3,,. 


and  the  eigenvectors,  orthonormallzed  with  respect  to  j|  JJ  and  the 

G 

corresponding  inner  product  are 


To  discuss  admissible  input  elements  in  this  case  we  let  Pk  ,  P  ^ 
be  non-negative  numbers,  k  =  1,  2,  3, ... ,  and  define 


H  {i«k>  *  3k  ,  M-  {-iuk}  =  P_k,  k  a  1,  2,  3,  ...  , 

GO 

M(Re(z)aa}-  U  ({ic^}  U  {-i«k»)  =  0. 


.  39) 
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Let  P(u),  -oo  <  <o  <  oo  t  be  defined  as  the  piecewise  linear  function 
such  that  In  the  interval  [  icj  ^ 1  “fc  +1  3 

\K+i-“)  + “  "k* 


3(«)  = 


(3.44! 


“k+1  "  “k 


SlftCe  S  V2i6  1/2 

r“k+l  1/2 a.  *  ,  _  KWk  +\+l"k+l  r  i 

J  «  '  P(  u)du  =  - 2 - l“k+l  "  “k-* 

wk 

—  4  (\  +\+i)  >  k  —  ®  , 


we  conclude  that  M-  is  a  Carle  son  measure  Just  in  case  there  is  a  constant 
C  such  that 

f  u^P(m)  du  =s  CJt-  or |  (3.45 

a 

whenever  0  <  a  <  t  ,  together  with  a  comparable  condition  involving  the 
and  negative  values  of  a>  .  But  (3, 45)  is  true  Just  in  case 

“lc  2  ’  k  «  1,  2,  3, . . . 

and  die  comparable  condition  for  negative  k  Is 

<*y/  s  C  ,  k  =  1,  2,  3, . . .  . 

Thus  for  the  Inhomogeneous  system 


(0-  *(0*(h)  ■ 

ft) 


the  Input  element 


with 
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Is  admissible,  from  this  criterion,  if 

k1/2(|Pk|2  +lP_k|2)  <  C  (3.46) 

for  some  fixed  positive  number  C  .  It  will  be  noted  that  this  is  (slightly) 
less  restrictive  than  the  requirement 


Bcample  5.  Negative  Results .  For  any  system  similar  to  the  one  in 

-“k 

condition  will  be  stronger  than  requiring  b  e  X  .  Hence  failure  of  the 
Carleson  measure  condition  cannot  be  used  to  show  that  an  element  b  is 
not  admissible ,  for  any  b  e  X  is  admissible  . 

To  illustrate  what  can  be  done  in  a  negative  direction,  we  return  to 

Example  1  with  0=0.  This  situation  has  been  studied,  using  a  different 
approach,  in  [13]  .  We  present  here  an  argument  more  in  the  spirit  of  the 

present  work.  As  shown  In  (3. 15), 


Example  4  but  with  |«k+1 


|  =  &(l/k'  +  e)  the  Carleson  measure 


\  2  v2  2 
\  =  -wk  s  "k 

and  (cf.  (3.16)  and,  w.l.o.  g.  taking  a  =  1 ) 


(3.47) 


v.% 


(3.4 


=  n2  kir  cos  ( kv )  -  ( -1 )  *v2  kir 


Since  *  j  bkj6  =  2k2  rr2  while  (k  +1  f  ir 6  -  k*  rr*  =  2k  ir*  +  ir*  . 

2 

it  is  not  hard  to  see  that  the  measure  y.  =  ,  p  ( -\k)  «  J  bfcJ  is 

not  a  Carleson  measure  in  this  case  .  As  we  have  remarked,  this  by  Itself 
is  not  enough  to  show  that  the  input  element  with  coefficients  (3. 48)  is  not 
admissible  .  To  show  this,  we  ask  the  reader  to  consider  the  function 

»|ir(2)=  (z  +l)’r 

analytic  in  the  complex  plane  minus  the  cut  consisting  of  {z  |  z  real,  z  =s  -1} 
If  r  >  ,  ipr  is  square  integrable  on  any  vertical  line  {z  |  Re  (z)  = 

£  ,  £  at  0}  with  uniformly  bounded  L  norm  and  i|/r(z)  is  bounded  for 

Re  (z)  0  .  It  follows  that  t|i  r(z)  is  the  Laplace  transform  of  a  function 


while  (k  +1  )2  ir 2  -  k2  t r2  =  2k ir2  +  ir2  , 


Re(z)  0 


fr=fr(t)  with  f r  €  I*  1 0,  ®  )  . 


Then 


'  k2  2 

e-K  ir  t  £^(t)  dt 


=  +r(k2.2,  . 

r  (kir  +l)r 

=  »(|k|1"2r),  k  -  -  . 

This  expression  is  not  square  summable  If  r  satisfies  the  inequalities 

l-2r  s-i  . 


so  we  require 


Vvf* 


Let  E  be  the  closed  subspace  spanned  by  the  functions  e 


-k2*2t 


In 


Lz[0,«)  and  let  ET  ,  T  >  0  ,  be  the  subspace  of  L^[0,T] 
consisting  of  restrictions  to  [0,  T]  of  functions  in  E.  If  T 
the  orthogonal  projection  of  f  onto  E  we  clearly  have 


is 


.2  2.  ~  »  .2  2. 

e-k  ir  t  f  (t)dt  =  /  e"k  v  1  f(t)dt. 

0  r  0  r 


jf 


It  is  shown  in  [7]  ,  [18]  that  the  natural  restriction  map  R  :  E  — >  ET 

is  onto,  (obviously)  bounded  and  (not  so  obviously)  boundedly  invertible 

2  2 

with  respect  to  the  induced  L  [o, «  )  ,  L  [0,T]  topologies  of  E,  ET  , 

,22.  1 

respectively.  Thus,  with  p^(t)=  e  , 


OB  ,2 

f  e-k  " 

0 


2t 


^{t)dt  =  (fr 


Pk>L2[0, 


00  ) 


=  R^Rjl  )  ,  =  «R~  V  f  ,  Rp.  )  , 

r  K  L2[0,»)  r  K  L£[0, 

<Pr(t )  dt 


T] 


-  /VkV* 


where 

9r=  (R_1)*£r  €  et  C  L2[0,T]  . 

It  follows  that  <pr  Is  an  element  of  L  [  0,  T  ]  such  that  the  numbers 

T  2  2 

(-l)k>JTkvf  e’k  *  Sr(t)dt,  k=l,2,3,... 

0  r 

are  not  square  summable  .  From  earlier  developments,  the  input  element  b 
with  coefficients  (3. 12)  corresponding  to  the  boundary  condition  (3. 10) , 
with  Ps  0  ,  a  =  1  : 


x(l,t)  *  u (t)  , 


is  not  an  admissible  input  element 


4.  A  Proof  of  Theorem  2. 4  .  It  is  clear  that  the  Car leson  measure  theorem 
in  ,  Theorem  2.4,  is  central  to  our  work  in  this  paper  .  This  result, 

in  one  form  or  another  has  been  known  for  somewhat  more  than  a  decade.  A 
proof  for  H^(D)  ,  where  D  is  the  unit  disc  in  the  complex  plane, 
appears  in  Duren  [  5]  .  A  proof  for  functions  in  is  given,  by 

Koosis  in  his  recent  book  [  12  ]  .  The  reader  is  also  referred  to  the  recent 
book  [8]  by  J.  Garnett .  Because  the  result  is  not  particularly  well  known 
outside  the  circle  of  mathematicians  working  in  Hp  theory  and  because  the 
results  are  rather  scattered  and  not  readily  available  in  precisely  the  form  we 
require,  we  offer  here  a  proof  of  Theorem  2. 4  which  is  a  direct  adaptation 
to  the  half  plane  of  the  result  for  the  unit  disc  appearing  in  Duren1  s  book  . 

The  proof  given  here  originally  formed  part  of  the  first  author4  s  doctoral 
dissertation  [9]  .  As  in  Duren*  s  work,  the  proof  makes  use  of  a  relatively 
simple  case  of  the  Marcinkiewicz  interpolation  theorem  ([23]  ,  Chapter  XI ) 
and ,  again  following  Duren,  we  do  not  quote  the  general  Marcinkiewicz 
theorem  but,  rather,  give  a  direct  proof  for  the  simple  special  case  required 


here  . 


We  begin  with  a  covering  lemma  of  "Vital!  type". 


Let  {I^|X.€A}=  J  be  a  family  of  intervals  In  R1  . 

Suppose  there  is  a  positive  number  K  such  that  for  any  finite  collection 

(1^  ,  1^  ,  •••  ,  1^  }  of  disjoint  intervals  in  J 
1  2  n  n 

E  IK  I  <  k. 

ks  1  k 


AJMJk* 


v;.\ 


■s 


Then  we  can  choose  a  sequence  {1^  |  ks  1,2,3,...}  of  disjoint  intervals 
from  3  with  the  property:  for  every  X  e  A  there  exists 
kc  {1,2,3,...}  such  that 


IXC^k 

is  the  interval  having  the  same  center  as 


but  five  times 


the  length  of  1^ 


Proof.  From  (4.1)  it  follows,  in  particular,  that  the  length,  |I^|,  of 

I ^  is  uniformly  bounded  (take  n  =  1,  X^  a  X  )  .  Define  the  sequence 

{I*  }  inductively  as  follows  .  Let  I.  be  such  that 
k  *1 


I'x.l  aLsu».  UJ  • 


For  k  =  2,  3, 4, 
and  such  that 


1  X  c  A 

let  I.  be  disjoint  from  I.  ,  i=  1, 2,  *..,n 
k  i 


|IX  |  *  |  sup{  1 IK |  |  X  e  A,  Ix  fl  \  ,  i  =  1,  2,«..,k-l>. 

lc  1 


Since  the  1%  are  disjoint  it  follows  from  (4.1)  that 


Um  |I\  I  =  0  • 
k  —  •  k 

Then  there  exists  k  such  that 
^  £  &  • 


Otherwise  (4.  2)  and  (4. 3)  could  not  both  be  true  .  Let  kQ  be  the 
smallest  integer  such  that  (4. 4)  is  true  .  Then 


0 


and,  together  with  the  foot  that 
IXC  Tk  ,  completing  the  proof 


this  implies  that 


We  subdivide  the  rest  of  the  proof  of  Theorem  2. 4  into  several  propositioi 
for  clarity  .  The  proof  is  given  for  the  half  plane  Re  (z)  >  0  ,  without 
loss  of  generality,  and  we  designate  H2  simply  by  H2  . 


Proposition  4. 2 .  Let  e  H2  and  let  ^(i  •  )  be  the  corresponding 

2 

boundary  function  in  L(-«,«).  For  z  =  a+ix,  a  >  0  ,  let 
I2  Ji®  the  interval 


and  let 


where  J 

z 


I2  =  [t-  a ,  T+  a] 

?(z)  s  sup  y|y  /|^0(it)|dt, 

is  the  set  of  all  finite  intervals  containing 

|*(z)|  *  £  fa). 


Then 


Proof.  From  the  Poisson  integral  formula  in  the  half  plane  we  have 


i  °(it)  dt_ 

W  +  (T-t)2 


so  that 


s  / 

LN=0  ?N 


,N+1. 


2  as  |t-xj  ^2 

a|j^0(it)jdt-| 


a|^Q(it)dt 
c2  +  (t-  t)2 


(4.5) 

(4.6) 

(4.7) 
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s  lT  y  r  l^Q<it)idt  +  li<o»t)|dt  -I 

*LM=0,t.T|s2N«a  4"o  a  J 

**[Z  -3TT  ?w**?wl 

Ln=o  2  J 

-  a.  *.> . 

fgopoaitlon  4.3.  Let  ip  e  L*(  -® ,  ®  )  and,  for  z  =  a  +  ir, 

a  >  0  ,  let  I,  be  given  by  (4. 5)  while  ( cf.  (4. 6)) 

$(z)=^sup  -jyy  ^  |4»(t)|dt  .  (4.8 

Let  |x  be  a  Carle  son  measure  and;  for  s  >  0  let  Eg  be  the  Borel 
measurable  subset  of  {z|Re(z)>0}  given  by 

Eg  =  {z  |  Re  (z)  >  0  ,  liT(z)  >  3}  , 

Then,  with  A  as.  bn  Definition  2.  3  , 

*  §  l+l.  ,  •  <*•’ 

L  (-»,  *) 

Proof.  Let  J  be  the  family  of  all  finite  intervals  in  R*  such  that 

JjJ  j f  |*(t)|dt  >  a  .  (4.1! 

If  Ij,  13,  •••,Ine  J  are  disjoint,  then  (4.10)  gives,  for  every  n, 

E  iys?E  /  i*(t)i«Jt 1  .  (4.« 

k  =  l  k*l  \  M-.-) 


Thus  J  satisfies  the  hypotheses  of  Lemma  4. 1  and  we  can  find  a  disjoint 
sequence  { ^  J  n  *  1,  2,  3, ...  }  c  J  such  that,  JR  having  the  same 


center  as  In  but  five  times  the  length,  each  Ic  j  is  contained 


some  Jn  . 


If  z  e  E.  then  I  c  I  for 


some  I  c  J  and  we  have,  for 


some  n , 


Then  clearly. 


lz  -  [t  -  t  +  cr]  c  J  . 
e.  n 


z  €  Sn  a  {  cr  +  iT|  0 


<  a  ^  ,  T  €  I  >  . 


This  being  true  for  all  z  e  E 

s  ' 


Es  c  U  -  S  . 
s  n  =  1  n 


Since  h  is  a  Carleson  measure  and  (4. 11)  holds, 

QB 

s  E  f*<Sn»=A  l 


n  =  1 


n=  1 


=  f  E  ig  =  H  i*i  ,  . 

n=  1  L  (_a,»  *) 

Proposition  4. 4  .  Let  c  Hg  with  boundary  function  (i  •  )  e  L^(-®, 

—  *(z)  be  defined  bjr  (4.6).  Then,  If  ,,  Is  a  Carleson  measure  on 

{  z  |  Re  (z)  >  0  }  , 

4(s)  >fiZ))Z  d|1  (Z>  *  10A/j*S0(lt>|2  *  .  (4. 


Proof.  For  each  r  >  o  let 


+r(t)  =  I 


»!0(lt)  If  |>i0(it)|  >  r 
0  otherwise  . 


z 

From  (1*  )  e  L  (-»,  «)  ,  we  conclude  that  the  support  of  4ir  is  a 
subset  of  (-«,  •)  of  finite  (Lebesgue)  measure  .  Then 

'I'r  €  ^  E  (  Tjr)  and  we  conclude,  since  t[<r  vanishes  outside 

£.r  *  that  e  <a)  ,  Moreover 

X  1+rl.l,  ,  dr  *  /  I  l>S0(lt)l<ltdr 

-  l>1o<lt>l 

=  /  /  =  /  l#so(lt>l  ^ 

-  ®  0  »  00 

■  K<‘*  >11*2,  • 

I^t  o(s)  a  p(Eg)  •  Then  we  can  see  that 

/  (?(z))2dp.  (z)  =  -f  s2da(s)  a  2/  sa(s)ds  . 

we  (z)  >  0  J0  0 

From  the  definition  (4.  8)  of  i|i  it  Is  clear  that  for  any  two  such 

functions,  tyj,  ^  ,  we  have 


Hence 


(+1  +  +2  >  C  z )  =s  ^(z)  +  $2(z)  . 

?(z)  *  (+r  +  (^0(i-)  -»l»r))(z) 

s  »Pr(z)  +  -  %)  (z) 

*  fr(«)  +  r 


■  ■'  t".  ,  "  ■.’V  V 


>1 

»sl 

I 

><l 

•  I 

T.I 


since  |^(it)-  4»r(t)J  is  either  equal  to  0  or  is  s  r  .  Le 

Fg  =  (z  |  ^(z)  >  s>  . 

Suppose  z  e  E2r .  Then  $  (z)  >  2r  and  (4. 15)  gives 

4>r  (z)  a  ?  (z)  -  r  >  r 
and  we  conclude  z  e  Ff  .  Thus 

E2r  C  Pr  • 

Hence,  from  (4.9)  of  Proposition  4.  3, 

M  V  *  MFr/2)  =  T^I+rU  , 

'  L  (  -  00  f  CO  ) 


so  that 


/  r«(r  )  dr  =  /  rjx(E  )  dr 
0  0  r 


15  5A  f  |{  «|i J  i  dr  =s  (using  (4.13)) 

0  r  L  (-«,  ®) 


*SAl^0(i-)l% 

L  (-<*>,  ob) 

Then  (4. 14)  gives  the  inequality  (4. 12)  . 

The  proof  of  Theorem  2. 4  is  completed  by  combining  (4.7)  of 
Proposition  4. 2  with  (4. 12)  above  to  give 

/  ,  .  I^(z)|2  dp(z)  as  /  (?(z))2d|i(z) 

Re(z)>0  v  Tte(z)>  0 


1000A 


/  |*0(it)|2dt 


»N  **•  i'V*  '*'•  - 


•  •••**«  »  *  ■  ‘  .V  ' n  \  \ 
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as  claimed  in  (2.  26) ,  except  for  the  trivial  detail  of  replacing  0Q(i  •  ) 

by  +a  -  I  ( a  +  i  *  )  . 
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